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Preface 
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existed where the non-rotating planet equations of motion 
were equivalent to the rotating planet equations of motion. 
This thesis is the result of this curiosity. 
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Rodney Bain, for his enthusiastic assistance, his patience, 
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Notation 


Roman Letter Symbols 


Acceleration (ft/s^) 
Ballistic coefficient 
Drag coefficient 


Lift coefficient 


Drag (lb) 


g - Acceleration of gravity (ft/s ; 

g - Gravitational acceleration at the planet’s 

* 2 

surface (ft/s ) 

h - Non-dimensional altitude 

I - Orbital inclination (deg, rad) 

L - Lift (lb) 

L/D - Lift-to-drag ratio 

M - Mach number 

m - Vehicle mass (lb ) 

v m 


Cosine of the flight path angle 
Radius from center of planet (ft) 
Planetary radius (ft) 


Aerodynamic reference area (ft 2 ) 

Time (s) 

Speed ratio, a modified Chapman variable 

Velocity (ft/s 2 ) 

Altitude (ft) 
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IX 
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•i 


Greek Letter Symbols 


Longitude of the ascending node (deg, rad) 
Inverse atmospheric scale height (1/ft) 
Flight path angle (deg, rad) 

Small parameter 
Small parameter 


Small parameter 


Longitude (deg, rad) 


3 2 

Planet graviational parameter (ft /s ) 


Subscripts 


Magnified non-dimensional altitude 


Density (lb m /ft; 


Bank angle (deg, rad) 

Latitude (deg, rad) 

Heading angle (deg, rad) 

Planet rotation rate (deg/s, rad/s) 
Argument of latitude at epoch (deg, rad) 
Angular velocity (deg/s, rad/s) 


At the surface of the planet 


Superscripts 


Inner 


Composite 


Outer 
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Abstract 


The assumption of a non-rotating planet, common in most 
analytical entry trajectory analyses, has been shown to 
produce significant errors in some solutions for the lifting 
atmospheric entry of Earth. This thesis presents an 
investigation of the validity of the non-rotating planet 
assumption for general three-dimensional Earth atmospheric 
entry. 

In this effort, the three-dimensional equations of 
motion for lifting atmospheric entry are expanded to include 
a rotating planet model. A strictly exponential atmosphere, 
rotating at the same rate as the planet, is assumed with 
density as a function of radial distance from the planet’s 
surface. Solutions are developed for the non-rotating Earth 
equations of motion and for one of the rotating Earth 
equations of motion using the method of matched asymptotic 
expansions. 

It is shown that the non-rotating Earth assumption 
produces incorrect entry trajectory results for entry 
orbital inclination angles between 0.5 and 75.0 degrees and 
vehicle speeds ranging from circular orbital velocities to 
low supersonic speeds. However, a variety of realistic 
trajectory states exist where some of the non-rotating Earth 














equations of motion are valid. Three of the non-rotating 
equations of motion are found to be valid for the same entry 
trajectory states. Other, independent trajectory states 
exist where a fourth non-rotating Earth equation of motion 
is valid. A fifth equation of motion is never valid for the 
ranges of orbital inclination angle and speeds investigated. 
Trends in the results of the trajectory states of validity 
are discussed and methods to estimate some of these states 
are presented. 









INVESTIGATION OF THE VALIDITY OF THE 
NON-ROTATING PLANET ASSUMPTION FOR 
THREE-DIMENSIONAL EARTH ATMOSPHERIC ENTRY 


I. Introduction 


Analytical studies often have the advantage of 
displaying valuable solution trends, giving more insight to 
the problem and its solutions than corresponding numerical 
work. Simple and accurate analytical methods to find 
solutions to the equations of motion for high speed flight 
vehicles supplement more complex and unwieldy numerical 
methods. In past analytical work on lifting atmospheric 
entry, the limiting assumptions of planar entry and a non- 
rotating planet were common. The equations of motion for 
planar atmospheric entry of a non-rotating planet are highly 
nonlinear; adding rotating planet and non-planar effects to 
the equations of motion make them even more difficult to 
solve analytically. Hence, the current literature has no 
investigations which analytically solves the equations of 
motion for three-dimensional, lifting atmospheric entry of a 
rotating planet. 

The Non-Rotatine Planet Assumption 


Although common in most analytical entry trajectory 
analyses, the assumption of a non-rotating planet model can 
produce significant errors in trajectory results. Since 
















most planets have a significant rotation rate, the rotating 
planet assumption will give more accurate entry vehicle 
performance results, especially for lifting vehicle range 
and time of flight calculations. Inherent in the concept of 
a rotating planet for atmospheric entry is the assumption 
that the planet’s atmosphere rotates with the planet at a 
constant rate. This assumption is more accurate than the 
non-rotating planet/atmosphere assumption but is still not 
ideal. The atmosphere of a rotating planet can be viewed as 
a boundary layer with rotation rates which change with 
altitude. Near the planet’s surface the atmosphere rotates 
at about the same rate as the planet. As altitude 
increases, the atmosphere rotates with a decreasing rate, 
and at high altitudes this rate is near zero. Hence, the 
true effect of a rotating atmosphere is therefore bounded on 
one end by the non-rotating planet solutions and on the 
other end by the rotating planet solutions. In this study 
is is assumed that the rotating planet solutions are ideally 
valid. 

The maximum effect of the rotating atmosphere (Chapman, 
1959:5) on a flight vehicle is easily calculated for a 
spherical planet. This maximum occurs at the equator for 
minimum altitude and is given by the ratio of the planet’s 
angular velocity to the circular orbital velocity at the 
planet’s surface. The planet’s rotational velocity on the 
equator is given by 










where u is the planet’s rotation rate and r is the radius 

of the planet. The circular orbital velocity at the surface 
of the planet is 

V = (g r ) 5 /2 

where g^ is the gravitational acceleration at the surface of 
the planet. The ratio is given by 



This ratio gives an indication of the possible error 
introduced to velocity calculations for a non-rotating 
planet model. For Earth, this ratio is approximately 0.06 . 
For Mars, Jupiter, Saturn, and Venus, this ratio is 
approximately 0.07, 0.30, 0.40, and 0.0, respectively (Vinh 
and others, 1980:3). Hence, for near-equatorial atmospheric 
entry, the maximum velocity error is negligible for Venus, 
significant for Earth and Mars, and very large for Jupiter 
and Saturn. The velocity error associated with the non¬ 
rotating planet assumption causes an even larger error in 
the calculation of aerodynamic forces. Since lift and drag 
are proportional to the square of velocity (as seen in 
Section II), the error introduced to the aerodynamic force 
calculations can be a maximum of about 0.14, 0.12, 0.60, and 
0.80, for Mars, Earth, Jupiter, and Saturn, respectively. 












A recent study of AOTV (Aero-assisted Orbital Transfer 
Vehicles) trajectories about Earth investigated possible 
trajectory simulation error due to the non-rotating planet 
assumption. In this study it was found that the non- 
rotating planet model caused velocity errors which gave 
dynamic pressure differences ranging up to 10 - 14% 

These differences may cause underprediction of the final 
altitude and overprediction of the attainable orbital 
inclination change in a non-rotating Earth analysis. It was 
concluded that rotating Earth effects must be included for 
realistic AOTV trajectory simulation (Ikawa, 1986:1,9). 

Another recent study (Miller, 1986:14) noted significant 
errors in values of range and time of flight for equilibrium 
glide entry trajectories when the Earth’s rotation was 
neglected. Differences in trajectory results between the 
rotating and non-rotating cases were found to be 
significant, especially for trajectories beginning at speeds 
greater than 15,000 ft/s. 

Vinh (Vinh and others, 1980:3) argues that inclusion of 
a rotating atmosphere into an analytical entry study would 
cause unwarranted overcomplication to the problem, possibly 
overshadowing the advantages of an analytical versus 
numerical analysis. However, for many atmospheric entry 
trajectories, such as multiple pass aerobraking, synergistic 
orbital plane change, and general high L/D vehicle 
trajectories with shallow entry flight paths, trends in the 
solutions caused by the rotating planet and its atmosphere 





















may be important. For these types of trajectories, the 
error in the calculation of the aerodynamic forces is more 
prominent due to the relatively large flight times within 
the sensible atmosphere. 


Problem 

Because of their complexity, the equations of motion for 
three-dimensional, lifting atmospheric entry of a rotating 
planet have not been analytically solved. However, the 
assumption of a non-rotating planet, common in most 
analytical entry analyses, has been shown to produce 

i 

significant errors in some solutions for the lifting ( 

atmospheric entry of Earth. An investigation of the general 

validity of the non-rotating planet assumption for three- 

dimensional Earth atmospheric entry is needed. In addition, 

the existence of trajectory states where the rotating planet 

terms in the equations of motion are negligible should be 

investigated. This would indicate the existence of 

trajectory states where existing solutions to the non- 

rotating equations of motion are valid for Earth lifting 

atmospheric entry. 

Scope 

In this effort, the three-dimensional exact equations of 
motion for lifting atmospheric entry are expanded to include 
a rotating planet model. The rotating planet terms in the 
equations of motion for Earth atmospheric entry are 
examined. Solutions are developed for the non-rotating 



















equations of motion and for one of the rotating Earth 
equations of motion. This was accomplished by treating 
atmospheric entry as a boundary layer problem, and applying 
method of directly matched asymptotic expansions. A variety 
of realistic Earth entry trajectory states are shown to 
exist where some of the non-rotating equations of motion are 
valid for a rotating Earth. This validity is coordinate 
dependent since sigularities exist in the equations of 
motion. Entry trajectory state examination is limited to 
orbital inclination angles between 0.5 and 75.0 degrees, 
where most Earth atmospheric entry occurs, and vehicle 
speeds ranging from circular orbital velocity to low 
supersonic speeds where terminal maneuvers, such as landing 
approaches, are usually initiated. 


Assumptions 

The planet is modelled in this analysis by a sphere 
having a central gravitational force field obeying the 
inverse square law. A strictly exponential atmosphere, 
rotating at the same rate as the planet, is assumed with 
density as a function of radial distance from the planet’s 
surface. In this effort, the only forces acting on the 
lifting vehicle are assumed to be gravity, lift, and drag; 
magnetic, solar wind, and other forces are considered 
negligible. The lifting entry vehicle is modelled as a 
point mass in a three degrees-of-freedom analysis. Constant 
lift-to-drag ratio is assumed along the trajectory and a 


















ballistic coefficient is specified for each flight 
vehicle/atmosphere under study. Angle of attack was not 
explicitly modelled. More detailed discussion of the 
approximations and assumptions will be presented in Section 
III. 

Approach 

In Section II, the equations of motion for three- 
dimensional lifting entry for a spherical, rotating planet 
are derived. In Section III these equations of motion are 
transformed into a form more convenient to examine and 
solve. The equations are also made non-dimensional and a 
coordinate transformation is undertaken. In Section IV the 
rotating planet terms in each of the five equations of 
motion are examined. These terms then are set equal to zero 
and checked for the existence of real solutions. It is 
shown that three of the equations of motion have identical 
solutions for these rotating terms. It is also shown that 
real solutions do not exist for the rotating terms in one of 
the equations of motion. The solution to the equation of 
motion containing these rotating terms is developed in 
Section V along with the solutions to the non-rotating Earth 
equations of motion. These solutions are derived from the 
rotating Earth equations of motion using the method of 
matched asymptotic expansions. In Section VI, the solutions 


to the rotating term equations in the other four equations 
of motion are examined in more detail. Trajectory states 











are presented where some of the non-rotating Earth equations 
of motion are independently valid for a rotating Earth. 
Methods are given to estimate solutions where four of the 
non-rotating equations of motion are valid for rotating 
Earth entry. Conclusions and recommendations for further 
study are presented in Section VII. 























In this section, the equations of motion are derived for 
three-dimensional lifting entry of a rotating planet. A 
spherical, rotating planet model is employed and it is 
assumed that the atmosphere rotates at the same rate as the 
planet with rotation rate, u. The lifting entry vehicle is 
modelled as a point mass in a three degrees-of-freedom 
analysis. Gravity, lift, and drag are assumed to be the 
only forces acting on the vehicle; magnetic, solar wind, and 
other forces are assumed to be negligable. Further 
discussion of assumptions and approximations is presented in 
Section III. 

Coordinate Systems 

Figure 1 defines the planet centered coordinate systems 
used in this analysis. The planet’s inertial reference 
frame has coordinates X, Y, and Z, with unit vectors I, J, 
and K, respectively. The X and Y axes lie in the planet’s 
equatorial plane and the planet rotates about the Z axis. 

The rotating planet-fixed coordinate system has axes X Q , Y Q , 

and Z Q , with unit vectors I Q , J Q , and K c , respectively. 
Another view of the planet centered coordinate systems is 
illustrated in Figure 2. a is a unit vector in the XY 


plane, pointing radially away from the planet. 6 is defined 







































Figure 2. Coordinate Systems Revisited 


as the longitude and is measured from the X Q axis which 

rotates about the K Q direction with rate u. f is latitude 
and is measured positive from the equator to the pole in the 
K q direction. 7 is the flight path angle and is measured 

"positive up" from the local horizontal to the velocity 
vector. The heading angle, is measured from a constant 
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systems. The velocity of the vehicle with respect to the 
inertial frame, X, Y, Z, is defined as the sum of the 
velocity of the vehicle with respect to the rotating frame 
and the cross product of the angular velocity of the 
rotating frame and the radius vector. 




IJK 


dr 

dr 


dt 

IJK dt 

+ fi x r 

lift 



( 2.10 


where r is the radius vector, extending from the planet’s 

center to the flight vehicle, ^ is the velocity vector, and 

^ is the angular velocity of the rotating frame. 

The inertial acceleration, a, is defined as the 
derivative of the inertial velocity. 
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d 2 ? 
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This can be simplified to the following 

















where 


IJK 


is the acceleration of the vehicle with respect 


to the inertial planet centered coordinate 
system. 



is the acceleration of the vehicle with respect 


to the radius vector, 
is the angular velocity 

a point mass vehicle in 
planet, is equal to a 


of the rotating frame. For 

flight over a rotating 
constant, the planet’s 


rotation rate, u 



is the Coriolis acceleration. 


a x 



is the Centripetal or Transport 
acceleration. 


To apply Eq (2.12) in the derivation of the equations of 
motion, transformations between the various coordinate 
systems are required. The coordinate transformation from 
the inertial system, XYZ, to the vehicle centered system, 
xyz, is not difficult but is prone to algebraic error 
because of the many intermediate steps and variables 
involved. For ease in derivation, this transformation of 
coordinate systems is divided into a number of simple angle 
rotations. 
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1st Rotation, Inertial to Rotating . The first rotation 
is the planet rotation about the Z axis. The rotation angle 
at time t is wt. From Figure 1 it can be seen that the 
following expressions apply: 


I Q = Icoswt + Jsinwt 
J = -Isinwt + Jcoswt 


K 0 = K 

In matrix form 


lo 

J 


coswt sinwt 0 

-sinwt coswt 0 


I 

J 

. K o. 


0 0 1 


K 


and 


I 


coswt -sinwt 0 


lo 

J 

= 

sinwt coswt 0 


lo 

K 


0 0 1 


K o 


( 2 . 1 ) 


( 2 . 2 ) 


2nd Rotation, Longitude . The second rotation (Figure 5) 
is a longitude angle rotation about the Z axis: 


a = I cos© + J sin© and j = -I^sin0 + J^cos0 

o o o o 


3rd Rotation, Latitude . The third rotation is a 
latitude angle rotation about the y axis. From Figure 6 


I = acos^ + K sin^ 


and 


Jc = -asin^ + K cos0 


The results from the longitude and latitude rotations are 
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combined and presented below: 
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and 


cos^cos# cos^sin# sirup 


?o 

-sin# cos# 0 



-cos#sin0 -sin#sin^ cos <p 


K o 



cos<pcos6 

-sin# 

sin^cos# 


♦ 

1 

-sin6cos<p 

cos# 

sin^sin# 


* 

3 

-sirup 

0 

cos^ 


k 


(2.3) 


(2.4) 


At this point enough information is known to make 

coordinate transformations between the XYZ inertial frame 

and the rotating xyz vehicle point mass frame. However, the 

transformation between the xyz and the X 0 Y 0 Z 0 frames is 

required. Figure 4 presented the vehicle centered 

coordinate systems that are used in the equations of motion. 

6 is a unit vector that is used to make the intermediate 

coordinate transformations easier to follow. 

4th Rotation, Heading Angle . The fourth rotation is a 
heading angle rotation about the x axis. From Figure 7 


B = jcosj} + ksin^fr and k 0 = -jsin^ + Jccos^ 

5th Rotation, Flight Path Angle . The fifth rotation is 
a flight path angle rotation about the z Q axis. It can be 

seen from Figure 8 that 
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The results from the heading angle and flight path angle 
rotations are combined and presented below: 
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Combining all five rotations gives the relations for a 
complete coordinate transformation between the x 0 y 0 z 0 and 

XYZ systems. 
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where c is a unit vector in the direction of the lift force. 
From Figure 6 it can be seen that aerodynamic, lift acts 

in the c direction, drag acts in the -j direction and the 

vehicle’s velocity vector points in the +j direction. 


Hence, lift is given by 

£< = Lc = (Lcosa)I 0 + (Lsina)k 0 

From Eqs (2.9) and (2.5) 


£ = Lcoscr ^cos'll - sin'ycos^j - sin^sin^Jcj 
+ Lsincr [-sin^j + cos *t] 


(2.13) 


£. = L(cosacos 7 ) J - L(cosasin^cos^ + sinasin\t)j 

+ L(sinacos^ - cosasin^sin^)lc (2.14) 

The drag components are found in terms of the vehicle 
centered unit vectors by application of Eq (2.6): 

S = -Dj q = -D(sin 7 )i - D(cos7cos^)j - D(cos 7 sin^)Jc (2.15) 

The vehicle’s local or relative velocity with respect to its 
own reference frame is 



V(sin 7 )i + V(cos7cos^)j + V(cos 7 sin^) lc 


(2.16) 


The planet’s rotation needs to be accounted for to 
obtain the vehicle’s inertial velocity from Eqn (2.16). A 
velocity component due to the planet’s rotation is added to 











the vehicle’s local velocity to form the inertial vehicle 
velocity. The planet’s rotation velocity component is a 
function of latitude; at the equator this velocity component 
has its maximum value and at the poles it is zero. 

Vehicle’s Inertial Velocity = P - 

IJK 

Velocity component due to planet rotation = ^ 

Vehicle’s local velocity = 

ijt 

$ ~~~ = $ + $ (2.17) 

IJK ijk r 

= V r j = wrcos^j (2.18) 

Combining Eqs (2.16), (2.17), and (2.18) gives an expression 

for the inertial velocity of the flight vehicle. 

^-= V(sin 7 )i + V(cos^cos^) + wrcos^) j 

IJK 

+ V(cos7sin^i) k (2.19) 

Another expression for inertial velocity can be derived 
from Eq (2.10) and compared to Eq (2.19) to produce three of 
the equations of motion. 

A dr dr * „ -» 

IJK dt IJK dt Ijfc I3R 

In this equation tl is the vector sum of two rotation 

rates between three coordinate systems, the XYZ, X Y Z , and 

J ’ 0 0 0 













xyz systems. Therefore, fl is the sum of the angular 
velocity of the rotating planet frame about the inertial 
frame and the angular velocity of the vehicle centered 
reference frame about the rotating planet frame. This can 
be expressed by 


a = a 


XYZ—»X Y Z 

0 0 0 


XYZ 

0 0 0 


►xyz 


From Figure 1 


- £ + i + $ - [w + H]k 0 - If 5 


and with Eqn (2.4) 


= [ w + <it] sin ^ - It 3 + [ u + ff] cos ^ 


( 2 . 20 ) 


Since 



ri 


dr 

dt 


m 


dr ? dl _ dr? 

dt 1 + r dt »•*«■ dt 1 


and 



T [ w + af ] cos ^ 3 + It 6 


( 2 . 21 ) 


Hence, the second equation for inertial velocity in terms of 
the vehicle centered reference frame coordinates is 




= If 1 + rcos^(w + |£)j + r|£fc 


( 2 . 22 ) 





















Three kinematic equations of motion are derived by 
equating the two inertial velocity expressions, Eqs (2.19) 
and (2.22), and comparing like terms. 


+ rcos <p(w + 4£) j + = V(sin 7 )i + 


+ V(cos 7 cos rf> + urcos^) ] + V(cos 7 sin^) k 
Kinematic Equations of Motion. 


i terms: 


St = Vsin ^ 


(2.23 


j terms: -r- 


VcosTfcos; 

rcos^ 


(2.24 


k tei 


Vcosnfsin! 


(2.25 


Derivation of the Force Equations of Motion 

To derive additional equations of motion, the inertial 
acceleration is calculated. From Eq (2.12) 


= a + 2fl x r 

ilk 


+ x r + ft x fixr 

ijk ljk ljk 


Each of the four terms in the above equation is derived 
separably for clarity. 


1st Acceleration Term. a 


From Eqs (2.21) and (2.23) 













Therefore a = si 

ijfc L 


in 7^ + Vcos 75?J i 


2nd Acceleration Term. X r 


From Eqs (2.20) 


j* d01 . ,* d <t> -> d0 ] ,r 

n = u + smpi - + u + cospk 


and from above 


-* dr» 

r *•»(• _ dt 1 
ijk 


dr? , r . ? 

= TT 1 = Vsin 7 i 


Taking twice the value of the cross product of these two 
equations gives the second inertial acceleration term. 

2?l x r = |2Vsin 7 cos^ (w + ^)]j + 2Vsin 7 ^jk (2 

ljk 


3rd Acceleration Term. ^ x r 


Taking the time derivative of Eq (2.20) gives 


* r del ,d^r d 7 e . ,? . r . dei . ,di d 7 <t>, 

= l w+ dtj^dt 1 + + r + dtj sin ^dt - 


■ [ w + Jt] sin 4t fc + TT 006 ^ + [" + 5t] cos ^H 

at 


Rewritten 



























d e 


Therefore, the third term of Eq (2.12) 


4th Acceleration Term. fl X fl X r 


Taking the cross product of Eqs (2.20) and (2.21) gives 


X ? = r [ w + ff] cos ^ + r ft* 

1 J K 


Taking another cross product produces the fourth inertial 


acceleration term. 


ft x (ft x r ) 

*56 


r d 01 2 2 . 

r [ w + dtj cos * 


d 2 ^ * 
+ r—l 


r ft[ W + ff] sin ^ 3 + r[u> + sin^cos^ k 


Total Inertial Acceleration . Combining Eqs (2.26), 


(2.27), (2.29), and (2.30) gives the total inertial 


acceleration in terms of the vehicle centered direction unit 


vectors. 


a-= Is 

IJK l 


. dv v di r Aey 2 . a 2 4 * 
in 7dt + Vcos 7dt - r r + dtj cos * - r 7T 1 

at/ . 






SKM 
















+ 2Vsin 7 cos^w + ^-J + r——cos^ - 2]u) + j^-Jrsi 

L dt 


r [w + sin^cos^ + r—^ + 2Vsin 7 ^ k 


(2.31) 


This equation can be simplified by substitution of the 
three kinematic equations of motion that were previously 
derived, Eqs (2.23), (2.24), and (2.25). Because of the 
size of this new equation, each direction component of the 
inertial acceleration, Eq (2.31), is examined separably 
below. 

Ith Inertial Acceleration Component . The Ith term of Eq 
(2.31) is 

. dV d7 T del 2 2 . d 7 <b 

ith = sm 7 ^ + Vcos 7j £ - r[w + cos 0 - r~f 

at . 

Substituting in for ^ and ^ with Eqs (2.24) and (2.25) 
Ith = sin 7 j£ + Vcos 7 jJ - rw 2 cos 2 0 - 2u;r [ Vc °^° s ^ ] cos 2 <p 


o 1/2 2 2 

2 . V cos nrcos 


,,2 2.2 
V cos nrsin 


r cos f 
This can be simplified to 


dV d'l V 2 2 

= sin 7 ^ + Vcos 7 ^f - —'cos ^ - 2u/Vcos 7 cos^cos0 


2 2 

u rcos <p 
















This equation simplifies to the following expression 


2 

j th - ^ C °cosd) S ^ TCOS( P + 2wVsin7Cos <p + 2^’sin7cos^cos if) 


2 

- 2 wVsin^cos 7 sin|t - 2— ’tan^cos 7 sin^cos^ 


The first term of this equation is 


d Ycoe'ycosit , ,dV . ,d 7 

dt " rcos 4> TCOS< f > = cos 7 c osV'dt " Vsin 7 cos^ 


,dtf VcoeTfcostf ,dr . ,d <t> 

- Vcos 7 s in ^-r 2/ cos *dt " rsin *dt 

r cos q> 


Simplifying and substituting in Eqs (2.23) and (2.25) gives 
3t VC rcosi S * rcOS< f > = cos 7 cos^ - Vsin 7 cos^ - Vcos 7 sin^ 

























V 2 V* 2 

- -■ sin^cos'ycos^ + -'cos ^cos^sinjttan^ 


Therefore, the entire jth term of Eq (2.31) is 
jth = cos7cos^ - VsinTfcos^^ - VcosTfsin^g^ 

+ 2wV^sin'Ycos0 - sin^cos7sin^J 


s^[si 


— * cos'jcosf |^sin7 - COS7SI 


in^tan^j 


(2.33) 


kth Inertial Acceleration Component. The kth term of Eq 


(2.31) is 


lcth = r ^ + 2 Vsin 7 ^^ + r sin^cos^ 

The first term in this fcth component requires some 
manipulation. Differentiating Eq (2.25) and multiplying by 
r gives 


d 2 d> d 
r—= r • — 


fcosTfsinii . ,dV ,, . . ,d7 

^-*■ = cos^smy^- - Vsin^siny^^ 


+ Vcos^cos^^ - sin7cos7sin^ 


(2.34) 


Substituting Eqs (2.34), (2.23), and (2.24) into the lcth 
term produces 


,V 2 


2 2 o 

V cos nrcos ifrsinij 


f- . 1 _V. ., I . , , V COS T cos 

kth = 2— • sin'jcos'jsxnf + w rsinpcosp + - rcosft — 

dV d'Y 

+ 2 wVcos 7 cos^sinrf + cos'isinit-rr - ’Vein'isiurb-rr 


















+ Vcos'ycos^Ji _ Bin'jcos'jsinip 


This equation can be simplified to the following expression 


\r 2 dV 

kth = -■ COS7(COS7COS jttan^ + ein^sini/i) + cos^sin^j^ 

- Vsin^sin^^ + Vcos'ycos^^ + 2u/Vcos7cos^sin0 


+ u 1 rsin^cos^ 


(2.35 


T otal Inertial Acceleration . Adding together the 
ith, j th, and kth acceleration terms given by Eqs (2.32), 
(2.33), and (2.35) gives an expression for the total 
inertial acceleration on the vehicle. 

-» dV d'Y V 2 2 

a-= sin'll + Vcos'f^f- - 2 uVcos' 7 cos^cos^> - — * cos 7 

- w 2 rcos ^ <p 1 + cos 7 cos^^- - Vsin 7 cos^^ - Vcos^sin^^^ 
+ 2Vw (cos<psin r ) - sin^cos^sin ifi) 


+ —•cos^cos^(sin^ - sin^cos^tan^) 3 


• j.dV v .dtf 

+ cos^sin^ + Vcos^cosy^-^ 


- Vsin7siny/^ + 2u/Vsin^cos7cos^ 


(2.36 


2 * V . 2 c 

+ u rsin^cos^ + —*C0S7(sin7sin^ + cos ^tan^cos7) k 






















Another equation for inertial acceleration can be 
derived by examining the forces on the vehicle. It is 
assumed that the only forces acting on the vehicle are 


gravity and aerodynamic lift and drag. 


Therefore 



and 


IJK 


d$ 

dt 


IJK 


= ;(£ - 6 ) 


+ e 


(2.37) 


(2.38) 


where the gravitational force, g, is a function of r and 
acts in the negative radial direction. 


g = ~g(r)i (2.39) 

Lift and drag are given by Eqs (2.14) and (2.15): 

t . = (Lcosctcos7) 1 - (Lcosasin 7 cos^ + Lsinasin^)j 
- (Lcosasin 7 sin^ - Lsinacos^)lc 

5 = (-Dsin 7 )i - (Dcos 7 cos^) j - (Dcos 7 sin^) 1c 

With substitution, a second equation is found for the total 
inertial acceleration of the vehicle: 


IJK 

[L D . ]t 

—cosacos 7 - ^sm 7 - gl 1 



1 

[L • _ ^ . L . 

D 



- ^cosasin7sin^ - ^sinacosjj + j^cos7sin^J & (2.40) 







By setting the two different equations for inertial 


acceleration equal and comparing like terms, three 
additional equations of motion can be derived. Equating Eqs 
(2.40) and (2.35) produces the following expressions: 


I Direction Term. 


dV 

sln 7dt 


d'V V 2 

Vcos^j-^ - 2 u/Vcos 7 cosyicos^ - —•cos 7 


2 2 

u rcos <p 


-(Lcosacos'j - Dsin 7 ) - g 


(2.41) 


3 Direction Term . 


cos 


dV . d 7 

S7~p£ - Vsinnf^f - 




n 7dt " Vcos 7 tan ^at 


2Vw 


+ (cos^sin7~sin$cos7siny) + 


V 2 

— * COS7 (sin7-sin^cos7tan^>) 


-1 

m 


(Lcosasin 7 


+ Lsinatan^ + Dcos 7 )j 


(2.42) 


Jc Direction Term. 


dV 

cos 7dt + 


wCosj dj£ 
tan^i dt 


Vsin 7^t + 


2wV si ^ c ° s 7 

tany 


+ 


w 


si 
r— 


n^cos^ 

sin^ 


V" 

+ — • COSTf 

r ' 


sin7 + 


cos^tan^cosnf 

tan^i 



(Lcosasin 7 


y sing 
^tan^ 


+ 


DCOS 7 ) 


(2.43) 


These three coupled equations can be reduced by some 
manipulation. Multiplying Eq (2.43) by -1 and adding the 
product to Eq (2.42) gives 


- Vcos7tan 



+ 


2Vw 

COBlfl 


(cos^sin7 - 


sin^cos7sin^) 












- ? -cos 7 [sinycos 7 tan^] - VffjJJ - 2^ V si g^° t 

_ w 2 r sin^cos^ _ V 2 . cos rcos^tan^cosj ] 
sint t> r ' L tanti J 


cosytan^cos' 

ta,nip 


L . , , L sina 

= - —sinatantt - — •--r 

m T m tan rp 


By combining terms and noting the trigonometric 


identities 


2 2 _ 1 2 2 

1 + tan ^ = (cos ip) and sin f + cos y = 1 

the above expression can be rewritten to form a new equation 
of motion: 


L. sing 
m cos 7 


+ 2Wui (cos0tan 7 sin^ - sin^) 


- lj r 


cos 7 


- —•cos 7 cos^tan^ 


(2.44) 


Substituting Eq (2.44) back into Eq (2.43) produces 


dV dT V 2 

0 = cos 7 ^ + 2Vu;cos0sin 7 cos^ - Vsin 7 j^ + -• cos 7 sin 7 

+ ^cosasin 7 + j^cos 7 + w 2 rsin^cos^sin^ 


(2.45) 


Multiplying Eq (2.41) by - and adding this 

product to Eq (2.45) gives another equation of motion: 


dtl 51 ^ + 


cos 2 Trl V 2 rcos 2 nr . 1 

sin 7 J r ' lsin 7 'J 


L F . cos 2 

+ —cosg sin 7 + —:- 

m l ' sin 7 


Ur 1 ] " + 2Vu/cos^cos^[sin 7 



















Rewritten, 


? S + w 2 rcosi sin^sintf + cos^^?^] 
sin^ J r L ~ ~ r sm7J 


= ^*cos7 + “cose - gcos^ + 2Vwcos0cos^ 
+ w 2 rcos^ |sin^sin^sin , y + cos^cos^j 


( 2 . 46 ) 


Substituting Eq ( 2 . 46 ) into Eq ( 2 . 42 ) produces the last 
equation of motion: 


dV 2 L 2 n \r * 

sin7j^ + —’cos 7 + — c.os'jcosa - geos 7 + 2 Vujcos')cos<pcosf 


u/ 2 rcos^cosip cos0cos7 - sin^sin^sin7j 


v 2 2 2 L 

- 2wVcos7cos^cos^ - —* cos 7 - w rcos ^ - - , cosacos7 


+ ^sin7 + g 


This can be rewritten 


dV D . 2 r 1 

= - — - gsin7 + u rcos 0 Icos0sin7 - sin^sin^cos7 1 ( 2 . 48 ) 


The Equations of Motion 

In summary, the following six equations of motion for 
three-dimensional atmospheric entry for a rotating planet 
have been derived and are listed below for convenience. 

^ - gsin7 + w 2 rcos^fcos^sin7 - sin^sin^cos7j ( 2 . 49 ) 













= — 'COS') + — • COSIT - geos'! + 2Vwcos^cos^ 
+ w ? rcos^i |sin^sin^sin')f + cos^cos'jfj 

vj* = ^ + 2Va/|cos^tan 7 sin^ - sin^>j 


u 7 r 


sind>cos(ftcosV* 
cos'y 


- -■ cos'ycos^tan^ 


dr 

dt 


- Vsin7 


d6 

dt 


dt 


V cos^cosip 
rcos <p 

Vcos^sin rj> 
r 


In the next section, these equations of mot 


ion are 


transformed into a form more convenient to analyze. 



















Ill. Approximations and Manipulation 
of the Equations of Motion 


The equations of motion for three-dimensional 
atmospheric entry for a rotating spherical planet were 
derived in Section II. In this section, assumptions and 
approximations used in this thesis are defined and 
discussed. In addition, the equations of motion are 
transformed into a form more convenient to examine and solve 
in later sections. The independent variable is changed from 
time to non-dimensional altitude, h, and the equations of 
motion are transformed into non-dimensional form by the 
introduction of non-dimensional variables. A coordinate 
system transformation is undertaken to utilize variables 
which are more convenient for atmospheric entry analysis. 

The equations of motion derived in Section II are 
given by 

dV D 2 

-rr = - — - gsin^ + u rcos^(cos^sin7 - sin^sin^cos7) (3.1) 


vdj V 2 
^dt = r cos 7 


+ u rcos <p ( 


ydj£ _ L. sing 
dt m COS7 


^cosa - gcos7 + 2Vu/cos^cosyi 
sin^sin^sin7 + cos^cos7) 

V 2 

-• cos7cos|ftan^ + 2Vu/(cos^tan7sin^ 


(3.2) 


sin^) 


0 

u rsinicostfcosf 




















dr ,, . 
dt = Vsin ^ 


d 6 _ VcoBTfcos{ 
dt rcosi 


d<t> _ Vcos^fsin; 
dt r 


(3.4) 


(3.5) 


(3.6) 


Note that the equations of motion for three-dimensional, 
non-rotating planetary entry can be derived from the 
rotating equations simply by setting the planet rotation 
rate, w, to zero. The equations of motion for the non¬ 
rotating planet assumption are therefore 


dV D 

dt = " m - esm7 


= £ 


>7 + -coscr - gcos 7 


2 

,,dt£ L sina V , 

y~rr = — ■-- - — • cosnrcosifltan d> 

dt m cosnr r i r r 


dr . 

dt = Vsin ? 

d0 _ Vcos7cos| 
dt rcos<p 

d d> Vcos 7 sint 


(3.7) 


(3.8) 


(3.9) 


(3.10) 


(3.11) 


(3.12) 


As expected, comparison of Eqs (3.1) - (3.6) and Eqs 
(3.7) - (3.12) shows that the equations of motion are 
significantly more complicated when the earth’s rotation is 
accounted for. 

























A few more variables require definition at this point. 


Let y be altitude and r + be average planetary equatorial 

radius. Since r is the radius measured from the center of 
the spherical planet to the flight vehicle, h is defined as 
the non-dimensional altitude and is given by 

h = (3-13) 

★ 


where 

r = r^+ y = r^+ hr^ = r^(l+h) (3.14) 

Therefore, = r + (3.15) 


, dt _ d r/dh _ r * 
dh dr/dt Vsin'y 

These relations will be used 
on the following pages. 


(3.16) 


in the approximations discussed 


Assumptions and Approximations 

The rotating planet assumption and the errors associated 
with the non-rotating approximation were discussed in 
Section I. Further assumptions and approximations are 
presented below. 

Spherical Planet Assumption . The approximation of an 
oblate planet by a sphere is very common in analytical 
flight mechanics analyses and is used here. It is a 
reasonable assumption for planets having small ellipticity 
such as Earth, Venus, and Mars but may not be as reasonable 
for planets such as Jupiter and Saturn which have relatively 















large equatorial bulges and ellipticities a magnitude 
greater than Earth’s (Chapman, 1958:2). For either case, 
error introduced by the spherical planet approximation is 
generally small for near equatorial entry trajectories. 

Spherical Atmosphere Assumption . The spherical planet 
approximation leads to the assumption that the atmosphere is 
spherically symmetric about the planet. In reality, planets 
are oblate spheroids causing their atmospheres to also have 
an oblate form. In addition, other significant deviations 
from the spherical model occur. A diurnal density bulge 
occurs over part of Earth’s sunlit side due to solar 
heating. Solar storms and fluctuations in a planet’s 
magnetic field can cause significant changes in density for 
a given altitude (Wiesel, 1986:66-69). However, since these 
effects generally occur at altitudes where aerodynamic 
forces are minimal (and inclusion of more sophisticated 
density models may introduce overcomplication of the entry 
problem) these effects are generally assumed to be 
negligable. The approximation of a spherical planetary 
atmosphere is perhaps the most limiting assumption to the 
atmospheric model. This approximation is better for the 
Terrestial planets, with their slow rotation rates, than for 
the large, outer planets (Duncan, 1962:276). 

Gravitational Model . The spherical planet approximation 
also leads to the assumption of an inverse square 
gravitational field. This is given by 














where is the gravitational acceleration at the planet 
surface. 

Atmospheric Density Model . As previously discussed, the 
planet’s atmosphere is assumed to be spherical and to rotate 
with the planet with a constant rotation rate, w . The 
planet atmosphere is mathematically modelled by an 
exponential atmosphere with the inverse atmospheric scale 
height, /J. This is a very common atmospheric model that has 
been successfully utilized in many studies on planetary 
atmospheric entry. Atmospheric density, p, is given by 
(Chapman, 1959:4) 

^ = -/Jdr (3.18) 

P 

It is sometimes assumed that the product of the inverse 
of the atmospheric scale height and the vehicle’s distance 
from the planet’s center is constant for a given atmosphere. 
With this model 

^ r or p = (r/r^) (3.19) 

where p is the density at the surface of the planet. 

The product fir has been approximated for many of the 
planet atmospheres. It has values on the order of 1000 for 
most of the planets; the mean value of /Jr is approximately 














900 and 350 for Earth and Mars, respectively (Vinh and 
others, 1980:5). 

For this study a strictly exponential atmospheric model 
is employed where the the inverse atmospheric scale height 
is assumed to be constant. This approximation allows the 
atmospheric density equation to be written in the form 

P = P* e ~^ y = P* e ~ h//e (3.20 


where e is a small number given by 



(3.21 


Note that is dependent on what planet is studied; for 
Earth, the average scale height is about 7.1 kilometers 
(Vinh and others, 1980:5) or about 23,300 feet. 

For planetary atmospheric entry, r is approximately 
equal to r^. This is an accurate approximation because the 

thickness of an atmosphere is generally very small compared 
to the planet’s radius. For example, the upper altitude 
limit of Earth’s sensible atmosphere is often taken to be 
350,000 feet. This value is only about 1.795 of the Earth’s 
radius. Approximating r by r^ leads to 


1 = _1_ 

P T * P T 

and hence the values given in the literature for mean 
planetary /Jr are considered equivalent to f)r . e is 


therefore a very small number (approximately equal to 1/900 


















for Earth) allowing it to be utilized as a small parameter 
in the asymptotic expansions of Section V. Earth 
atmospheric density values were calculated with this model 
and plotted in Figure 10 with values obtained from the 1976 
U.S. Standard Atmosphere (NOAA, 1976:Table IV). Comparison 
of these results shows that the accuracy of the exponential 
density model is reasonable for Earth with a constant value 
for /Jr of 900. 

Aerodynamic Forces and the Ballistic Coefficient . Lift 
and drag accelerations are given by the following familiar 
expressions: 


L P SV * C L , D ^ SV2C D 

m “ 2m and m ~ 2m 


(3.22) 


where m is the mass, S is the aerodynamic reference area, V 


is the velocity, C^ is the lift coefficient, and C^ is the 
coefficient of drag for the flight vehicle. 

The non-dimensional ballistic coefficient, B, is defined 
to help place the equations of motion into a form easier to 
solve. It is given by 


P~ SC T 


(3.23) 


and is specified for each flight vehicle under 
consideration. The ballistic coefficient is a function of 
the vehicle’s physical characteristics and the planet 
atmosphere and is considered to be a constant (Busemann and 
others, 1976:18). 
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cos'y (1+h) J 




Differentiating Eq (3.18) with respect to h gives 


du _ 2V(1+h) cos 2 7 . dV 


dh 


g r 
° * * 


dh 


2 • 2 2 
2 V (1+h) cosTfsinTr . d^ V cos i 

g r dh g r 

* 4c 4e 


(3.26) 


Substituting Eqs (3.20), (3.23), and (3.25) into Eq (3.22) 
gives a non-dimensional equation for the aerodynamic drag 
acceleration: 


— - v 2 —Qf ^/e - fl rv 5 » ^/e 


= V SC n e 
m 2m D 


= pBV e 


or 

D ug^Bexp(-h/e) 

m e( 1 +h)cos 2 7 

Likewise 

L = D.^L = C L u g»Bexp(-h/e) 
m m Cjj Cp e ( 1 +h) cos 2 7 


(3.27) 


(3.28) 


Changing the Independent Variable 

The equations of motion for a rotating planet are now 
rewritten using the expressions derived above. To change 
the independent variable from time, t, to non-dimensional 
altitude, h, the equations of motion are multiplied by Eq 
(3.16). Direct substitution into the dV/dt equation gives 
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O. 0. 1-3 Cl CL 

cr < cr tr < 


1 12 


1 /l 


ug^r^BexpC-h/e) 


e( 1 +h)cos^sii^ 


( 1 +h) 


ug r 


g r COS 7 


( 1 +h)' 


( 1 +h) 


ug r 

* 


2 2 ( 1 +h) 

+ w r —t — L - 

* sin^ 


( 1 +h) 


ug r 


1 / 2 


•cos0cos7(cos^sin^ - sin0sin^cos7) 


is equation can be reduced to 


Bexp(-h/e) 


ecos7sin7 


ug r 

* 


( 1 +h) 


1 / 2 


COS 7 


K* r * 


1 / 2 


2 COS $ 
+ U ~-*- ‘ 


tan7 


(l+h) 3 ^ 3 ! 1 /2 




u(l+h) J J 

(cos^isin 7 - sin0sin^cos7) (3.29) 


Direct substitution into the d 7 /dt equation gives 


a_.. 


C L u g»Bex P ( _ h/e) cosa r^(1+h)cos 2 7 


cos 
















€003731117 


(l+h)sin7 tan7 ug 


2 2 

u r^(l+h) sin^cos^cosf 

•(cos0tan7sin^-sin0) - ---- (3.31) 

u ®* tan / 


The d6/dt equation becomes 


dh (1+h)cos^tan7 

and the dip / dt equation becomes 

dip _ sin rp 

dh (l+h)tan 7 

The dr/dt equation is incorporated into the other 


(3.32) 


(3.33) 


equations, reducing the number of equations of motion from 
six to five. The du/dh equation also incorporates the dV/dh 
equation. Combining Eqns (3.19) and (3.20) allows the du/dh 
equation to be rewritten as 
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d 


2 


Substituting in for the dV/dh and d'y/dh terms using Eqs 
(3.29) and (3.30) eliminates V and gives the du/dh equation 
as a function of u. 


-tJexp ( -l 


(1 + h) ^ r g ecos 7 sm 7 (1+h) 


u e* r * 


cos^ 


.u (l+h)‘ 


tan7 ug 


(l + h) 3 r 3 1 1 / 2 


• (cos0sin7~sin$siny/cos7) 


2utan7 (]+h)tan7 + CL. esi 


L Bcos a -h/e 


(1+h)utan7 


+ 2u> 


(1 +h) r 1 1 


tan7 


r (1+h) . 

2 * COSffiCOS'Y \ 

+ w ----‘-(sinmsinwsinTf+cos'ycosm) 

ug tan7 v r r 1 ' r 


This rather formidable equation reduces to the following 


du -u 
dh - (1+h) 


2uBe \ L cosa 

- — ; + __- 

e sm 7 CL 0037 


- 4 u 


(1+h)ur 


* ._ j. _ l o. .2 /1 . i_ \ 2 * 


* costpcosip + 2w (1+h) — cos^/cos#' 


cosmsinTf sinmsintficos' 


tan7 


tan7 


- sin^sin^sin7 - cos^cos7 










a o- ^ n q-icl v c '- J ^ a a 

tr c W ►—* crUs « » rt> CD sr c 


The du/dh equation can be further reduced to 


-u 2uBe h/£ [, C L , 

ITThy - esirry [ 1 + tan 7 cosaj - 4w 


(l + h)ur 1 1/2 


• costpcosip - 2 u 


( 1 +h; r , , . 

2 * cosfflsinfflsint 


tan 7 


(3.34) 


Another variable is introduced to simplify the equations 
usemann and others, 1976:19). 


t q = cos'y h< 


dq . dnr 

dh = - Sin ^dh 


(3.35) 


ing these relationships, the d'y/dh equation can be 


written 


= _3_[s 2 _ jl _ ft. 

(l+h)lu *J C D 


c L Be - h /€ rr,(l + h)l ’/ 2 

~ -cosa - 2qwcos$cosii - 

C D e U Z* 


u ? r^q 3 ( 1 +h) 2 "^^(sin^sin^tan^ + cos^) 


(3.36) 


The equations of motion for three-dimensional, rotating 
anetary entry have now been transformed from Eqs (3.1) - 
. 6 ) to the following: 


2uBe 


-h/e 


( 1 +h) esin 7 


1 + *tan 7 cosa - 4w 


( 1 +h)ur 


•cosAcostf - 2u 


2 < 1+h > r * 


>sin®smt 


(3.37) 


cosa - 2 qu/cos 0 cosd» 




















(sin^sin^tan7 + cos^) 


(3.38) 


- w 2 r q 3 (l*h) 2 -£2^ 

V / Ug^ 


djt _ ^L . Be ^' //£ sinq 
dh - Cp ecos^sin^ 


cosTrcos^tan# ^ 
(1+h)sin^ 


2u 

tan7 


(l + h)r 1 1/2 


u 8. 


• (cos^itan^sin^-sin^) - 


2 2 

w r^(l+h) sin^cos^cos^ 
ug tan 7 


(3.39) 


dO _ _ cos^ 

dh (1+h)cos^tan7 


(3.40) 


dip _ si n^ 
dh (l+h)tan 7 


(3.41) 


The Classical Orbit. Variables 

Up to this point, the equations of motion have been 
presented as functions of latitude, longitude, heading 
angle, and other variables. This form has been often 
utilized for atmospheric trajectory simulation by numerical 
integration. For ease in studying and in deriving solutions 
in later sections, these equations are placed in terms of 
the classical orbital elements, ft, I, and a. For a non¬ 
rotating, spherical planet, these variables are constants of 
motion for non-atmospheric flight. This characteristic 
greatly simplifies the solution derivation for the non¬ 
rotating planet case. I is defined as the orbital 
inclination angle, fl is the longitude of the ascending node, 
and a is the argument of latitude at epoch. Basic spherical 
trigonometric relations are found in Appendix A and these 
relations are applied to transform the variables 6, $, and ^ 
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V. V, 



















in terms of the orbital elements, a, fl, and I. Figures 11 
and 12 show the geometry of the two sets of variables. 

The following relations, derived in Appendix A, relate 
6, <p , and ip, and a , 0, and I. 


sirup 

= sinlsina 

(3.4 

cosl 

= cosipcosip 

(3.4 

si nip 

tan# 
tan a 

(3.4 

sin (0 

tanji 
” tanl 

(3.4 

sin^ 

= sinlcos(0-fi) 

(3.4 

COSQ 

= cos#cos(0-fl) 

(3.4 


Differentiating Eq (3.43) gives: 
sinl'dl = cos^sin^'d^ + sin^cos^’d^ 
Therefore 

dl _ cos^sin^ . dip sin^cos^ . dip 
dh sinl dh sinl dh 


dl _ cos#s imp . dip sin^costftsinV' 
dh - sinl dh sinltan 7 (1+h) 


Substituting for d^/dh, Eq (3.39), gives 


-h/e 


dl _ _L, cosasinqBe ' _ cosqcosnfcost/’tan^ 

dh Cp ecos 7 sin 7 (l+h)sin 7 


+ 2 u 


r (l+h)l l/2 


u e. 


^^^(cos^tan7sin^ - sin^) 


(3.4 


2 

2 (1+h) . ^ ^ , sintficoslsintf , 0 ,, 

- u r •— k —i— z — cosasinmcosmcostp + -— utt- •—t (3.4 

* ug tan 7 r r r cos0(1+h)tan 7 sinl v 























igure 11. The Osculating and Equatorial Planes 
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The last term cancels with the second term and with 


application of some of the spherical trigonometric 


equations, Eqs (3.42) - (3.47), the dl/dh equation is 


transformed from a function of 6, <p , and ip , to a function of 


a , fi, and I. 


dl L cosasinaBe 


-h/e fr (l+h)l ,/2 


dh (L. ecos 7 sin 7 


+ 2 ui 


cosasinl , 
tan7 


(3.50) 


• (cosatanor - sina) - w 2 r *— ^ —(sinlcoslcosasina) 
v 1 ' * ug^tan7 v ' 


To find an expression for da/dh in terms of the desired 


variables, Eq (3.42) is differentiated. 


d <f> _ coslsina . dl sinlcosq _ da 
dh ~ cos <p dh cos^ dh 


Substituting for d<p/dh with Eq (3.33) produces 


da _ _ 1 _ _ tana . dl 

dh ~ (l+h)tan 7 tanl dh 


(3.51) 


da _ _ 1 _ _ ^L . sinasinaBe 

dh (l+h)tan 7 ecos 7 sin 7 tanI 


- 2 w 


r^(l+h) 


sinacosl 


tan7 


(cosatan 7 - sina) 


2 ( 1 +h) . 2 2 -r 

+ u r — J —sin a cos I 
* ug tan 7 


(3.52) 


To find an expression for dfl/dh in terms of the desired 


variables, Eq (3.46) is differentiated. 


sin^ = sinlcos( 6 - 0 ) 





















By substituting in Eqs (3.32) and (3.45), this 


expression can be rewritten as 


d# _ cosltan# . dfi _ cosltan# coslcosq . dl 

dh cos ip dh (1+h)tan^cos# cos^cos^ dh 


(3.53) 


Rewriting Eq (3.48) gives another expression for dip/ dh: 


d# _ sinl . dl _ sin#cos# 
dh - cos^sin ip dh (1+n) tan^ 


(3.54) 


Equating Eq (3.53) and Eq (3.54) gives dfl/dh as a function 
of dl/dh. 

dfi _ cos ip cosltan# _ sin#cos#cos# 

dh cosltan# (1+h)tan^cos# (l+h)tan , y 



dl 

r sinl 

coslcosal 


+ dh 

Lcos#sin# 

cos#cos#J 


dfi _ _ cos# __ cosltanc 

dh cos!tan#tan#(1+h)tan# cos# 


coslsin^ 

2 * 
cos # 


dl . cos ip r< 
dh cosIsin#L< 


s#cosaJ 


11 T • 1 VWOU/N-WDW 

dh coslsin^Lcosa r 

This simplifies to the following 

dfi _ tana .dl 
dh sinl dh 


(3.55) 


Eq (3.55) can be "de-simplified" by substituting Eq (3.50) 
into this equation. 


, n C T n -h/e . q Tr (l+h)l l/2 

dfi _ L . Be ' sinasma 2w * v ' 

dh ~ Cp esin , )(cos7sinI tan7 ug^ 
















• [^sinacosatan^ - sin aj - ug tan 7 —'coslsin a 


(3.56) 


The du/dh equation, Eq (3.34), can easily be transformed 
to the desired orbital variables by application of the 
spherical trigonometric relations. 


u 

2uB 

-We. 

(l+h) 

esin 7 



ur *( 1 + h) 

1 / 2 

4wcosI 


— 


g * 



jsatan'yj 


2 2 

2 (l+h) cosasin Isina 

j r ' -»- L --- 

* g tan 7 


(3.57) 


The last equation of motion to be transformed is the 
dq/dh equation, Eq (3.38). The w 2 term of Eq (3.38) is 

- J 2 r^q 3 ^cos^sin^sin^tan 7 + cos 2 <p j 

This term can be rewritten as 

2 3 (l+h) 2 . . 2 Tj . 2 T • 2 1 

- w r q cosasinasin Itan 7 + 1 - sm Ism a 

* ug* ' J 

The dq/dh equation transformed into the orbital elements is 
dq q fq 2 ,1 C L Be _h ^ E „ [ r »( 1+h )1' ^ 

dh = TuhTlu " J j “ —i— cosa ~ 2q “ cosI [ U g. j * 

- w 2 r^q 3 + ^ sinacosatan 7 sin 2 1 + 1 - sin 2 Isin 2 aj (3.58) 


In summary, the equations of motion for three- 
dimensional rotating planetary entry have been derived for 
the independent variable, non-dimensional altitude, and some 













V 


convenient dependent variables including orbital 
inclination, longitude of the ascending node, and argument 
of latitude at epoch. These equations are given below. 

du - _ _H_ 2Bu . -h/ef, Si. cosataml 

dh " (1+h) esin 7 e [ C D cosatan 7j 


- 4wcosI 


ur^(1+h) 


2 r . (1+h) 2 cososin 2 Isinq f3 
* g tan 7 


da_a_[a 2 _ il _ 

dh " (1+h) Lu X J c D 


C L Be h/e „ T [ r *( 1 

C n e ug 


r (l+h)l ,/2 


2 3 (1 + 

- w r a -*■- 

* ug 


dl _ L . cosasingBe 
dh - CL. ecos 7 sin 7 


-h/e 


+ 2ui 


r (1+h) 1 1/2 


cosasinl 

tan 7 


(cosatan 7 - sina) 


.2 _ . (1+h) 


- u r •—^—r——(sinlcoslcosasina) 
* ug tan 7 v 


C T n -h/e . . „ f r (1 + h) 1 ^ 2 

L . Be ' sinasino 2u; * 

Cjj esin 7 cos 7 sinI tan 7 ug^ 


sinacosatan 7 - sin 2 aj - 


w 2 r^(1+h) 2 

——t -* cos Is in a 

ug tan 7 


i C, . . n -h/e 

_1_ _ L . sinasingBe 

(l+h)tan 7 Cp. ecos 7 sin 7 tanI 


r J 1+h ) 1 1 /2 sinacosl , . . . 

i - —r-(cosatan 7 -sina) 


XVvV.vW? 


[sinacosatan 7 sin 2 I + 1 - sin 2 Isin 2 aj (3. 


(3.6 


(3.6 









in acos I 


(3.63) 


+ u r 


* ug tan 7 


Singularities exist for these equations of motion for 
flight path angle values of 0.0 and 90.0 degrees and for an 
orbital inclination angle of 0.0 degrees. Special 
consideration must be taken when evaluating Eqs (3.59) - 




















The rotating Earth terms in each of the equations of 
motion are examined in this section. This examination is 
conducted for orbital inclination angles ranging from 0.5 to 
75.0 degrees and vehicle speeds ranging from circular 
orbital velocity to low supersonic speeds, where terminal 
maneuvers such as landing approaches are usually initiated. 
The rotating Earth terms are set equal to zero for each 
equation of motion and are then checked for the existance of 
real solutions. The existance of real solutions for any of 
these expressions indicates trajectory states (specific 
combinations of values of u, h, q, I, 0, and a along an 
entry trajectory) exist where that particular non-rotating 
Earth equation of motion is valid. In a later section, 
these trajectory states will be examined in more detail. 

The non-existance of real solutions to the rotating Earth 
terms in any one of the equations of motion indicates the 
corresponding non-rotating equation of motion will be 
invalid for any Earth entry trajectory. A solution that 
accounts for the Earth’s rotation is needed for any non¬ 
rotating Earth equation of motion that is invalid for the 
full range of inclination angle and speeds investigated. 
These solutions are developed in Section V by application of 
the method of matched asymptotic expansions. 




















The Second Small Parameter 


To begin, the equations of motion are placed in a more 
convenient form to work with. In Section III, £ was shown 
to be approximately equal to the ratio of the atmospheric 
scale height and the radius of the planet. A convenient 
second small parameter is introduced to replace the w terms 
in the equations of motion. This second small parameter is 
defined as the square of the ratio of the planet’s 
rotational velocity at the equator and the prograde, 
equatorial circular orbital velocity at the surface. The 
orginal small parameter, e, is re-labelled e^ and the second 

small parameter is labelled e. 






- € 2 (1 + h) (tan^sinacosasin I + 1 - sin Isin a) (4.3) 


C T D -h/e 

L Be l sinasma 


da 1 _L. Be 

dh (1+h) tan 7 e ^ 


cos'ysin'ytanl 


1/2 

n \ Q+h)l sinacosl/ , . % 

- 2 e A - *- -t -(cosatani - sina) 

.2 u tan 7 v ' ' 


+ ( h) ;. sin!acos!l 

2 utan7 


(4.4) 


= gL. Be~ h/e i singsina + 2 f ii±hll 1 X !5isa (cOBat - sina ) 

dh Cp e^ cos'jsin'jsmi L 2 u tan 7 x ' 


2 

x 4 ,v 2 .sin acosl 
— e„(l+h) 7 

2 utan^ 


, T C T D -h/e 
dl L Be ' l sxnacosa 


(4.5) 


dh c D Cj 


cos^sin^ 


L 2 u 


l/2 

+h) 1 ' . cosasinl , 

u J tan^ 


x , . \ „ m 2 cosasinacoslsinl 

• (cosatan'y - sina) - C 2 (1+h) --- - 


(4.6) 


It was previously noted that the latter three of the 


above equations are coupled 


dfi _ tan g . dj[ 
dh sinl dh 


da _ _1_ _ tana . dl 

dh ~~ (l+h)tan 7 tanl dh 


The Rotating Earth Terms 


(4.7) 


(4.8) 


The equations of motion for lifting atmospheric entry 
for a rotating Earth can be written as the sum of the 
equations for a non-rotating Earth and of the terms that 




























Nonrotate 


(1+h) Lu A J C 


D e i 


• cosct 


= - COSl 


1 Rotat 


(4.10) 


2 Q 2 2 2 

- e (1 + h) (tan-ysinacosasin I + 1 - sin Isin a) 


Argument of Latitude at Epoch . The rotating Earth 
equation of motion for argument of latitude at epoch, given 
by Eq (4.4), can be rewritten as 


= Nonrotate + Rotate 


where 


da = 1_L.Be 

dh Nonrotate (l + h)t.an 7 “ Cp £ 


C T n -h/e 

L Be ' l sinasina 


cos7sin7tanI 


I Rotat' 


1/2 

«r (l+h)l ' sinacosl, , . x 

= - 2[£ 2 - i -^] tan 7 “ ■»■“) 


* c (l^h ) ? - sin? ° cos21 
2 utan7 


(4.11) 


Longitude of the Ascending Node . The rotating Earth 
equation of motion for longitude of the ascending node, 
given by Eq (4.5), can be rewritten as 


= Nonrotate + Rotate 


where 





















dn 

dh 


Nonrotate 


L . Be ' l singsing 
Cp cos7sin , ysinI 


dfi 

dh 


Rotate 



. sing 
tan 7 


(cosgtan 7 


sing) 


,, . IX 2 .sin 2 gcosl 
£ ( 1 +n ) t 

2 utan7 


(4.12) 


Inclination Angle Equation . The rotating Earth equation 
of motion for orbital inclination angle, given by Eq (4.6), 
can be rewritten as 

= Nonrotate + Rotate 


where 


dl 

dh 


C L Be~ h/£ l .si_ 

Nonrotate C D £ t cos 7 si n7 


nacosg 


dl 

dh 


- 2 

[e -^11 

Rotate 

.2 u 


tan7 


(cosgtan 7 - sing) 


, n 2 cosgsingcoslsml 

e (1 + h)- t- 

2 utan7 


(4.13) 


Investigation of Real Solutions to the Rotate Terms 

Eqs (4.9) - (4.13) are the Rotate equations, those parts 


of the complete equations of motion which account for the 
rotating Earth. In the following pages, the five Rotate 
equations are each set equal to zero and then examined for 
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1 


From Eq (4.9) 


Rotate 


- 4^e 2 u(l+h)j cosl 


2 e 2 ( 1 +h) 


2 

2 cosgsin Isina 
tan 7 


Setting this equation equal to zero and solving for one of 
the five variables gives 


7 = tan 

Rotate = 0 


2 

-1 -Ccososinqsin I 
2cosI 


(4.1 


where c=p^%J ,/2 


(4.1 


Substitution of a few realistic combinations of u, h, 
and a indicate real solutions exist for Eq (4.14). 

Rotate Term Solutions for the Inclination Ancle 


Equation . From Eq (4.13) 


3K __= z [ e 2 X1 ? 1 ] '“t”'” 1 !' 05111 "’ “ sin “> 


Rotate 


- e 2 ( 1 +h) 


2 . cosasi nocosIsini 
utan^ 


Setting this equation equal to zero and substituting in Eq 
(4.15) for C gives 


O * T S X nO V . . -y -r 

zsinicosa cosa - . - = t -sinacosasinlcosl 

tan7J tan7 


(4.16) 


One set of trivial solutions to this equation is in the form 


sinlcosa = 0 for any C and 7 























k». 


The w term contribution to the dl/dh equation of motion is 
therefore equal to zero when sinl and/or cosa = 0. However 
due to the fact the equations of motion are singular for 
sinl = 0 as discussed in Section III, the trivial solution 
is actually 

cos a = 0 (a = . . . -3ir/2, -n/2, r/2, 3tt/2, ...) 

Nontrivial solutions to Eq (4.16) can be found by solving 
for one of the five variables. 


7 = tan 

Rotate = 0 


l sing + Ccoslsina 
2 cosg 


(4.17 


Substitution of a few realistic combinations of u, h, I, an 
a indicate real solutions exist for Eq (4.17). 

Rotate Term Solutions for the Argument of Latitude at 


Epoch . From Eq (4.11) 


I Rotate 


l/2 

2 ( 1 +h) ' ' . sing cosl 

^12 u tan^ 


(cosgtan'j - sing) 


2 9 

,, L \2 sin g cos I 

+ e ( 1 +h)- 7- 

2 utan^y 

Setting this equation equal to zero and substituting in Eq 
(4.15) for C gives 


o • t sing] 

2 smgcosl cosg - ^ an ~j = 


tan7J tan^ 


C .2 2 T 

--‘sin geos 1 


(4.18) 


One set of trivial solutions to this equation is in the form 


sinacosl = 0 for any C and 7 





























The u term contribution to the da/dh equation of motion is 
hence equal to zero when sina and/or cosl are zero. 

However, due to the range of inclination angle selected for 
study, cosl is never equal to zero and the trivial solution 


sina = 0 (a = . . . - 2 ?r, -ir, 0 , ir , 2ir , •••) 

Nontrivial solutions to Eq (4.18) can be found by solving 
for one of the five variables in this equation. 


7 = tan 

Rotate = 0 


-1 Ccoslsinq + sing 
2 cosa 


(4.19) 


Substitution of a few realistic combinations of u, h, I, and 
a indicate real solutions exist for Eq (4.19). 

Rotate Term Solutio ns for the Longitude of the 
Ascend i ng Node . From Eq (4.12) 


= 2 


1 / 2 


Rotate 


( 1 +h)!' sina, ^ N 

[ e 2 ^rj t^,( cosat * n 'i ~ sinQ > 


e 2 (l + h) 


■ 2 T 
2 . sin a cosl 

utan7 


Setting this equation equal to zero and substituting in Eq 
(4.15) for C gives 


0 . sina n ■ 2 cosl 

2sina cosa - t - = Csin a -——- 

tan7_ tan7 


(4.20) 


One set of trivial solutions to this equation is in the form 
sina = 0 for any C, I, and 7 


















Hence, the u term contribution to the da/dh equation of 
motion is equal to zero for a = ... -2ir, -* , 0, if, 2 k, 


Nontrivial solutions to Eq (4.20) are in the form 



-l 

Ccoslsina + sina 

7 

= tan 

Rotate = 0 

2cosa 


(4.21) 


Substitution of a few realistic combinations of u, h, I, and 
a indicate real solutions exist for Eq (4.21). 

Coupling of Three Rotate Term Solutions . The non¬ 
trivial solutions for the rotating term expressions are 
identical for the da/dh, df]/dh, and dl/dh equations, Eqs 
(4.17), (4.19), and (4.21). This can also be seen by noting 

that the dl/dh, dfl/dh, and da/dh equations are coupled by 
Eqs (4.7) and (4.8). 


dfi _ tana.dl 
dh sinl dh 


(4.7) 


da _ _1_ _ tana . dl 

dh (l+h)tan 7 tanl dh 


(4.8) 


These coupling relations can be written in terms of the 
Nonrotate and Rotate expressions: 


dfi 

dh Nonrotate 
+ Rotate 


tana.dl 
sinl dh 


+ 

Nonrotate 


tana . dl 
sinl dh 


Rotate 


(4.22) 


da 

dh Nonrotate 
+ Rotate 


1 

(1+h)tan^ 


tana,dl 
tanl dh 


Nonrotate 


tana .dl 
tanl dh 


Rotate 


(4.23) 











Since the Rotate terms contain u or eby definition, the 

first term of Eq (4.23) must be a Nonrotate expression. 
Hence, 


da 

dh 

In 


Rotate 

addition 


tang . dl 
tanl dh 


Rotate 


(4.24) 


dn 

dh 


Rotate 


tang . dl 
sinl dh 


Rotate 


(4.25) 


From Eqs (4.7) and (4.8) it can be seen that the dfl/dh 
and dg/dh equations differ significantly. However, 
comparing like terms indicates 


da 

dh 


Rotate 


cos 



Rotate 


(4.26) 


These results simplify the search for the trajectory 
states where the non-rotating Earth equations of motion are 
valid. Instead of detailed examination of the non-trivial 
solutions of all three of the da/dh, dQ/dh, and dl/dh Rotate 
equations, examination of only one of them is required. The 
trivial solutions for rotating term expressions for the 
dfl/dh and da/dh Rotate equations are the same, sina = 0, and 
the trivial solution for the rotating term expression for 


the dl/dh Rotate equation is cosa = 0. 

Rotate Term Solutions for the Flight Path Angle 













dh 


Rotate 




COS1 


o 

2 Q • 2 2 2 

- e 2 ( 1 +h) (tan 7 sinacosasin I + 1 - sin Isin a) 


Setting this equation equal to zero gives 


3 2 2 2 1 

0 = 2qA^cosl +q tan 7 sinacosasin I + 1 - sin Isin al (4.27 


where 


A = 
1 


T 1 /2 


e 2 (1+h)' 


(4.28 


Because singularities exist for I or 7 = 0.0 or 90.0 

degrees, q = 0 can not be considered as a trivial 

solution. To examine possible solutions to Eq (4.27), some 

additional variables, A., are defined for convenience. 

’ x 


let 

A 2 

1 . 2 T . 2 

= 1 - sin Isin a 

(4.29 


A 3 

. 2 

= sinacosasin I 

(4.30 


A 

A 

= 2 A^ cosl 

(4.31 


Using these variables, dividing by q, and remembering 
q = COS 7 , equation Eq (4.27) can be rewritten 

A + q 2 A ta.n'j + q 2 A =0 

*1 u Z 

2 

A + A cos 7 = - A sin 7 cos 7 (4.32 

H Z o 


Squaring this expression 

























2 


A 2 cos 4 7 + A + 2A A cos 2 7 - A 2 cos 2 7 + A 2 cos 4 7 = 0 
2 * 4 4 2 ' 3 '3 ' 


(A 2 + A 2 )i / 2 + (2A A - A 2 )i/ + A 2 = 0 

2 3 42 3 4 


(4.33 


, 2 2 

where 1 / = q - cos 7 


(4.34 


Eq (4.33) can be easily solved in the form 


Ax/ + Biv + C = 0 


where A = A ? 2 + A 3 2 , B = 2A^ A ? - A^ , C = A^ 2 , and 


-B * [b 2 - 4Acl 
" = -2A- 

The Rotate term solution for the flight path angle equation 
is hence given by 


-2A A + A 2 * [a. 


- 4A A A / - 4A. A 


At this point 1 / has been placed in terms of the 
variables A^ , A^ , A^ , and A^ . Looking at each of these 

expressions individually aids in the determination of 
general trends of the values of u . 

A Equation. It can be easily seen that A has a 

maximum value of 1.0 when sina = 0 , for any value of sinl 

within the inclination range of interest, 0.5 I 75.0 

degrees. A minimum of 0.0670 occurs for A when I = 75.0 















degrees and a = 90.0 degrees or 90.0 ± 180.0 degrees. 
Thus, A 2 is never negative. 


A Equation. A has a maximum value of 0.466 for 

—3- - 3 

I = 75.0 degrees and a = 45.0 degrees. A minimum value of 


-0.466 occurs when I = 75.0 degrees and a = -45.0 degrees 


A^ Equation . Since A^ is always positive, it can 


be easily seen A has a maximum value of 34.0 for 

4 


I = 0.5 degrees, u = 1.0, and h = 0.0 . A has a minimum of 

4 


0.0 for u = 0.0 . Minimum values of A for non-zero u occur 

4 


for I = 75 degrees. For a prograde orbit A is always 

4 


positive 


Solutions . For real solutions to exist in Eq 
(4.35), and hence for the rotating terms in the dq/dh 
equation to ever have a zero contribution, the expression 
within the square root must be positive. Factoring A 3 out 

of the root expression in Eq (4.35) leaves the following 
condition for the existence of real roots 


A 2 - 4[A A + A 2 1 > 0 

3 L 2 4 4 J 

where A ? and A^ are always positive and A^ 


(4.36) 


can be positive 


or negative over the range of inclination angle examined. 
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Eq (4.36) indicates real solutions to the dq/dh Rotate 
equation do not exist for realistic entry trajectories 
except for values of u corresponding to low speed, subsonic 
and near subsonic flight. Therefore, no real roots to Eq 
(4.35) exist for orbital inclination angles between 0.5 and 
75.0 degrees and vehicle speeds ranging from circular 
orbital velocities to low supersonic speeds. This result 
indicates the non-rotating Earth dq/dh equation of motion is 
invalid for Earth atmospheric entry for the ranges of 
orbital inclination angle and velocity investigated. 
Therefore, a first order solution (Section V) to the 
complete dq/dh equation is required for Earth atmospheric 
entry analysis. This solution is developed in the next 
section using the method of matched asymptotic expansions. 
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Solutions to the Equations of Motion 


Matched Asymptotic ExDansions 


In Section IV, the terms in the five equations of motion 


that account for Earth’s rotation were examined. It was 


determined that trajectory states exist for four of the five 


equations of motion where the rotating Earth terms give a 


zero contribution. For the fifth equation of motion, the 


dq/dh equation, it was shown that the rotating earth terms 


always have a non-zero contribution for realistic lifting 


entry trajectories with entry inclination angles between 0.5 


and 75.0 degrees and speeds ranging from circular orbital 


velocities to low supersonic speeds. Therefore, to 


adequately model Earth atmospheric entry, a solution to the 


dq/dh equation is required that includes the rotating Earth 


effects. This solution, along with the non-rotating Earth 


solutions to the other four equations of motion, is 


developed in this section, by application of the method of 


matched asymptotic expansions. 


Combining the Small Parameters 


In Section III, e was defined to be the ratio of the 


planetary atmosphere scale height and the planet’s radius. 


A second small parameter was introduced in Section IV to 


non-dimensionalize the rotating planet terms in the 


equations of motion. This second small parameter was 


defined as the square of the quotient of the planet’s 


rotational velocity at the equator and the prograde, 




























equatorial circular orbital velocity at the mean planetary 
radius. 


e. = a-nd e„ = 

l p r 2 


The general equations of motion for three-dimensional, 
rotating planetary entry were given in Section IV in the 


(l + h) 
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dh Cp e ) cos7sin7 
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dl _ . Be singcosa + ^ [ e (l+h) 1 1 cosasina . 

dh _ Cp e ] cos 7 sin 7 ^L £ 2 u J tan 7 

• (cosatan 7 -sina) - ( 1 +h)- utan7 - (5-6) 

Values for the two small parameters must be calculated 
and substituted into Eqs (5.1) - (5.6) for each planet 
studied. For the Earth, the value of the first small 
parameter is approximately 1/900 and the value of the second 
small parameter is approximately 1/289 . To simplify the 
following analysis, it is noted that the first small 
parameter is approximately equal to one third of the second 
small parameter. Hence, a new small parameter, e, is 
defined for Earth entry analysis as 

e = = e , ,/2= <v 3)1/2; e i = g2 and s = 3e2 (5 ’ 7) 


The Equations of Motion for Earth Atmospheric Entry . To 
obtain the equations governing Earth atmospheric entry, the 
first and second small parameters are replaced in the 
equations of motion by e and a constant. The resulting 
equations are used in the derivation of solutions for Earth 
entry: 
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It was previously seen that the latter three of the 


above equations are coupled: 

dn _ tang .dl 
dh sinl dh 


da _ _ 1 _ _ tana . dl 

dh (l+h)tan 7 tanl dh 


(5.12 


(5.14 


Eqs (5.8) - (5.12) describe the three-dimensional, 
rotating planetary entry of Earth. However, the generi< 


IVOvCV'IV?/'v.v.\-v.v v/.-'v. -■'•--.v'v 
















ij 





application of these equations to some of the other planets 
is not lost. As an example, consider the problem of Mars 
atmospheric entry. To make Eqs (5.1) - (5.6) applicable for 
Mars entry analysis, the first and second small parameters 
are calculated using Mars planetary constants. A new small 
parameter, e, differing from Earth’s by only a constant is 
then defined. For Mars, the first small parameter has a 
value of approximately 1/350 . The second small parameter 
is calculated from IAU defined constants to be approximately 
1/218. 

11 2 r * _ 1 
£ 1 " /Jr ~ 350 e 2 " W g^ “ 218 

The first small parameter is approximately equal to 
three fifths of the second small parameter. Therefore, a 
unique e is defined for Mars entry analysis: 

e = 1^7 = (e,) l/2 = (3e, /5) 


or 


e 2 and e 

2 



(5.15) 


To adjust the equations of motion for Martian 
atmospheric entry, the new first and second small 
parameters, as a function of e and a constant, are 
substituted into Eqs (5.1) - (5.6). Hence, the equations of 
motion for Mars entry will differ from those for Earth entry 
analysis by only a constant in each expression containing e . 
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The Method of Matched Asymptotic Expansions 


While many analytical methods have been applied to the 
problem of planetary atmospheric entry, one method remains 
relatively unexploited. This is the method of matched 
asymptotic expansions. From an entry vehicle’s viewpoint, 
planet’s atmosphere forms a boundary layer of density in 
space. Aerodynamic forces on the vehicle change from 
insignificant to dominating for a relatively small change i 
altitude. A variety of methods, including composite 
expansions, multiple scales, and matched asymptotic 
expansions, have been utilized to solve boundary value 
problems. However, the method of matched asymptotic 
expansions is more "versitle" and "effective" than these 
other methods for both linear and nonlinear problems 
composed of partial or ordinary differential equations 
(Nayfeh, 1985:257,258). Past papers describe successful 
applications of matched asymptotic expansions to solve 
limited flight mechanics problems involving lifting 
atmospheric entry. (Busemann and others:1976, Shi and 
Pottsepp:1969, Shi:1971, Shi and others:1971, Willes and 
others:1967). Most of these efforts used two-dimensional 
equations of motion and assumed a spherical, non-rotating 
Earth model except Busemann (Busemann and others, 1976), 
where the three-dimensional equations of motion were 
employed. In the following pages, the solution for the 
dq/dh equation of motion for a rotating Earth is derived. 

To obtain this solution, the solutions for the non-rotating 














equations of motion are first derived from the rotating 
Earth equations of motion. 

For a boundary layer problem, the method of matched 
asymptotic expansions gives two or more solutions, each 
valid in specific regions of the domain. These solutions 
have some overlap and can therefore be matched. The 
matching conditions allow for the construction of a 
composite solution which is valid over the entire domain. 

In this analysis, the entry problem is treated as a boundary 
layer problem with a solution developed away from the 
boundary layer, valid for exo-atmospheric flight, and a 
solution developed at the boundary layer, valid for 
atmospheric flight. These two solutions to the equations of 
motion can be independently applied in their respective 
domains. However, it is more convenient to have a solution 
valid for the full range of entry conditions. Hence, the 
inner and outer solutions are matched and a composite 
solution valid in both regimes is constructed. 


Outer Expansions 

The solution developed for the exo-atmospheric portion 


of the domain is called the outer solution. This solution 


is developed from asymptotic expansions of the equations of 
motion using the small parameter, e. The outer solution 
variables are denoted by the superscript "o". Straight¬ 
forward expansions of the outer solution variables are made 
and are as follows: 
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(5.16) 


n + e‘n. + e 2 n„ + 


a + £ 1 a + e 2 a + 

0 1 2 


nr + e 1 <> + e 2 nf + 

'0 '1 '2 


The solutions for lifting atmospheric entry of a non- 


rotating planet are of order £ . Adding the rotating planet 


model to the atmospheric entry problem causes additional 


terms of order £ and £ 2 to appear in the expansions of the 


equations of motion. These additional terms are 


relatively complex but require the solutions to be carried 


out to order £ in order to account for Coriolis 


acceleration on the flight vehicle. To also model the 


usually insignificant Centripetal or Transport 


acceleration, solutions would have to be carried out to 


order £ . The solutions developed in this study are to 


order e° for Eqs (5.8), (5.10), (5.11), and (5.12) and to 


order £ 1 for Eq (5.9), the dq/dh equation of motion. 


Solutions developed to order £ J act as a correction to the 














zero order solutions, a correction which accounts for the 


significant effects of a rotating planetary atmosphere. 

In the following pages, the equations of motion are 
expanded for small e using the straight-forward expansions 
given by Eq (5.16). Many of the common expressions that 
contain e in the equations of motion are expanded in more 
detail in Appendix B. 

The du/dh Outer Expansion . The outer expansion of the 
du/dh equation is given by 
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The dl/dh Outer Expansion . The outer expansion of the 
dl/dh equation is as follows: 
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The dfl/dh Outer Expansion . It. has been noted that, the 
dO/dh equation is a function of the dl/dh equation as well 
as a function of a and I. This dependence exists for the 
outer expansions of these equations as well. From Eq (5.13) 
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The da/dh Outer Expansion . It has been noted that the 
da/dh equation is a function of the dl/dh equation as well 
as a function of a, h, 7 , and I. This dependence exists for 
the outer expansions of these equations as well. From 



















Eq (5.14) 
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The outer expansion of the da/dh equation is 
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The dq/dh Outer Expansion . The outer expansion of the 
dq/dh equation is as follows: 
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Solutions to this set of differential equations are 
derived in Appendix C and are given below. Eq (5.30) is the 
outer zero order e solution to the dq/dh equation of motion. 
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e 1 Terms . The complete set of e 1 term outer expansion 
differential equations, derived above, are 
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A complete solution to the dq/dh equation requires that 
the solutions to both Eqs (5.18) and (5.26) be found. 
Solutions to these differential equations are derived in 
Appendix C and are given below. Eqs (5.32) and (5.33) are 
the solutions to the order e terms in the outer solutions of 
the du/dh and dq/dh equations, respectively. 
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where C , C , and C are functions of C and C , and C 

13 ’ 14 ’ 16 1 2 ’ 20 

is a constant of integration. Expressions for these 

constants are given in Appendix C. 

Order e Outer Solution for the dq/dh Equation of Motion . 
The order e outer solution for the cosine of the flight path 
angle is given by the substitution of Eqs (5.30) and (5.33) 
into Eq (5.16). 
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Inner Expansions 

The solution developed for the atmospheric portion of 
the domain is called the inner solution. This solution is 
developed from asymptotic expansions of the equations of 
motion using the small parameter, e. The outer solution 
variables are denoted by the superscript "i". Straight¬ 
forward expansions of the inner solution variables are made 
and are given below. 
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To derive the inner expansions, a new independent 
variable is required so as to force the equations of motion 
to focus on the boundary layer region. £, a magnified 
version of h, is the inner expansion independent variable. 
Hence, £ is defined as the magnified non-dimensional 
altitude and it is given by 


£ = h/e 2 or h = e 2 £ 


(5.36 


Therefore ^ = e 2 and by the chain rule, 
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The equations of motion for Earth atmospheric entry are 
now transformed from functions of h, represented by Eqs 
(5.8) - (5.12), to functions of the magnified non- 
dimensional altitude, £. 
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The latter three of the above equations are coupled 
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In the following pages, these equations of motion are 
expanded for small e using the straight-forward expansions 
of Eq (5.35). Many of the common expressions that contain e 
are expanded in detail in Appendix B. 

Th e du/d£ Equation Expansio n. The inner expansion of 
the du/dh equation is given by 
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'0 '0 

The dl/d£ Equation Expansion . The inner expansion of 
the dl/d£ equation is 

dl dl C L 
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This equation reduces to 
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The d0/d£ Equation Expansion . It has been noted that 
the dfl/d£ equation is a function of the dl/d£ equation as 
well as a function of a and I. This dependence exists for 
the inner expansions of these equations as well. 
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The inner expansion of the d0/d£ equation is 
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With substitutions, the dfl/d£ equation expansion becomes 
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The da/d£ Equation Expansion . It has been noted that 
the da/d£ equation is a function of the dl/d£ equation as 
well as a function of a, 7 , and I. This dependence 
exists for the outer expansions of these equations as well 
From Eq (5.40) 
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The inner expansion of the da/d£ equation is 
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With substitution, the following expression is obtained: 
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The dq/df Equation Expansion. The outer expansion of 


the dq/d£ equation is 
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Inner Expansion Solutions 

e° Terms . The complete set of e° term inner expansion 
differential equations, derived above, are as follows: 
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Solutions to this set of differential equations are 
derived in Appendix D and are given below. Eq (5.55) is the 
inner zero order e solution to the dq/dh equation of motion. 
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e 1 Terms . The complete set of e 1 term inner expansion 
differential equations, derived above, are as follows: 
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Order e Inner Solution for the dq/dh Equation of Motion. 


The order e inner solution to the cosine of the flight path 
angle is given by the substitution of Eqs (5.55) and (5.60) 
into Eq (5.35). 

q 1 = cJq + eq i + 0(e 2 ) (5.35) 

C T 

q 1 = p^-Be ^cosa + K + eK (5.61) 

°D 1 7 

At this point, the inner and outer solutions to the 

equations of motion have been derived to order e° . These 
are the solutions to the atmospheric entry for a non¬ 
rotating Earth. In addition, the inner and outer solutions 
to the dq/dh equation of motion have been derived to order 
e. These latter two solutions include the Coriolis 
acceleration on the flight vehicle and therefore account for 
the significant effects of a rotating Earth. In order to 
create composite solutions valid for both the inner and 
outer regions, matching is performed. 

Matching Zero Order e Solutions 

One of the fundamental rules of the technique of matched 
asymptotic expansions is that the outer expansion of the 
inner expansion solution is equal to the inner expansion of 
the outer expansion solution (Nayfeh, 1981:277,278). This 
condition allows matching of the solutions, reducing the 
number of unknowns in the solution equations. On the 
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following pages the inner expansions are taken of the e 


outer expansion solutions. Next, the outer expansions are 


taken of the e° inner expansion solutions. Finally, the 


derived expansions of expansions are matched and the inner 


solution constants of integration are expressed in terms of 


the outer solution constants of integration. 


Inner Expansions of Order £ Outer Expansion Solutions. 


From Eq (5.29), the one term outer expansion for u is 


(u Q ) = C/(l+h) 


2 2 

Noting that £ = h/e or h = e £ and rewriting this 


equation in terms of the inner variable, f, gii 


(u 0 )° = C/(l + e 2 0 


Expanding for small e gives the inner expansion of the 


outer expansion of the zero order e term for u. 


K>°r- 


(5.62) 


From Eq (5.30), the one term outer expansion for q is 


(q o )° = C /[2(1 + h) - C C 2 (l + h) 2 ] 


Rewriting this equation in terms of the inner variable 
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Expanding for small e gives the inner expansion of the outer 
expansion of the zero order e term for q. 
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1 / 2 


(5.63) 


From Eq (5.28), the one term outer expansion for H is 



Rewriting this equation in terms of the inner variable 
simply gives 


= c , 


Expanding for small e gives the inner expansion of the outer 
expansion of the zero order e term for Q. 


K>°f - 


(5.64) 


From Eq (5.27), the one term outer expansion for I is 


(I ) = C 

v 0 J b 


Rewriting this equation in terms of the inner variable 
simply gives 
















Expanding for small e gives the inner expansion of the outer 


expansion of the zero order e term for I. 


[‘V 0 ] 1 - ° 5 


(5.65) 


From Eq (5.31), the one term outer expansion for a is 
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Rewriting this equation in terms of the inner variable: 
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Expanding for small e gives the inner expansion of the outer 
expansion of the zero order e term for a. 
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(5.66) 


Outer Expansions of the e° Inner Expansion Solutions 


From Eq (5.55), the one term inner expansion for q is 


(q Q ) i = 0^' Be ^cosa + K ] 


2 2 • • 

Noting that £ = h/e or h = e £ and rewriting this 

equation in terms of the outer variable, h, gives 




cos? + K 











Expanding for small e gives the outer expansion of the inner 
expansion of the zero order e term for q. 


and hence ) *j = 
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From Eq (5.56), the one term inner expansion for u is 
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exp 
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Rewritten in terms of the outer variable 
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Expanding for small e gives the outer expansion of the inner 
expansion solution of the zero order e term for u. 
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From Eq (5.59), the one term inner expansion for I is 
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Rewritten in terms of the outer variable 
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Expanding for small e gives the outer expansion of the inner 
expansion of the zero order e term for q. 
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From Eq (5.57), the one term inner expansion for a is 

sinK. 

(aj 1 = sin 


sin(I Q ) 1 


Rewritten in terms of the outer variable: 


sinK 


(a l 1 = sin 1 

° sin(I ) X J 


Expanding for small e gives the outer expansion of the inner 
expansion of the zero order e term for a. 


(5.71) 
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From Eq (5.58), the one term inner expansion for fi is 
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Rewritten in terms of the outer variable 
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(n A ) 1 = -cos 


• cos (a ) 1 - 


Expanding for small e gives the outer expansion of the inner 


(5.72) 


expansion of the zero order e term for Q. 

. , o _ [ . , o 

(n o ) 1 j = -cos 1 cos (Oq) 1 ] /cosK 4 + K 5 

f il ° 

where [ (a Q ) I was given by Eq (5.66) 


6° Term Matching . As previously discussed, the outer 
expansion of the inner expansion solution is equal to the 
inner expansion of the outer expansion solution. 


Mathematically, 


[w l ]° - [«°f 

This is sometimes expressed (Lagerstrom, 1972:90) as 

lim , lim. x = lim. lim . x 
outL in J ml out J 


(5.73) 


where x=x + ex + e x + 

0 1 2 


The expressions derived for the outer expansion of the inner 
expansion solution and for the inner expansion of the outer 
expansion solution can be matched for q, u, a, I, and Q. 

Matching the Speed Ratio Equations. From Eq (5.73) 


K>f - k>°f 
















Therefore, from Eqs (5.62) and (5.69) 
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Matching the Flight Path Angle Equations 




Setting Eq (5.63) equal to Eq (5.57) gives 
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Matching the Inclination Angle Equations 
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Equating Eqs (5.65) and (5.70) gives 
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Matching the Argument of Latitude at Epoch 
Equations . 


K>r ■ K>°r 
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Equating Eqs (5.66) and (5.71) gives 
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Matching the Longitude of the Ascending Node 


Equations 
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Setting Eq (5.64) equal to Eq (5.72) 
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The inner expansion solutions, valid in the domain of h 


near the boundary layer where aerodynamic forces are 


dominate, have been matched with the outer expansion 


solutions, valid in the domain of h away from the boundary 


layer. The results of the matchings are used to reduce the 


number of unknowns in the solution equations. The constants 


of the outer expansions are determined first, often from 
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initial conditions, in a planetary entry problem. Hence 
is desirable to express the constants of the inner 
expansions, K, in terms of the constants of the outer 
expansions, C. From Eq (5.75) 

K , - [V < 2 - °, C 2 >] ,/! < 5 

From Eqs (5.74) and (5.75) 
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From Eq (5.77) 
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From Eq (5.78) 
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From Eq (5.76) 
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With Eqs (5.79) - (5.83) the inner solutions can be 


written in terms of the constants, C. 


Zero Order e Solutions to the Equations of Motion 

A' solution good for the entire range of h is desired for 
the convenience of not having to apply the two sets of 
solutions to one entry problem. The composite solution, 
valid everywhere in the domain, combines the inner and outer 
solutions that were valid only in certain overlapping 
regions of the non-dimensional altitude domain. This 
composite solution is given in the following equation by the 
method of matched asymptotic expansions. 

(y)°] = y° + y 1 - [(y) 1 ] (5.84) 

The composite solution is therefore given by the sum of 
the inner and outer solutions and the difference between 
this result and any expression common to both the inner and 
outer solutions. In addition, the inner expansion of the 
composite solution is equal to the inner solution, and the 
outer expansion of the composite solution is equal to the 
outer solution. Mathematically, 
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Speed Ratio Composite Solution . 
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Combining Eqs (5.29), (5.56), and (5.62) gives the composite 
solution for the speed ratio. 
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Flight Path Angle Composite Solution, 


(5.86) 


q° - (q)° + (q ) 1 ~ [ Cq)°] 1 = (q)° + (q ) 1 - [ (q) £ ] ° 

Combining Eqs (5.30), (5.55), and (5.67) gives a composite 

solution for the cosine of the flight path angle. 


q c = Cj /[ 2 ( 1 +h) - C^U+h) 2 ] 


C L „ -h/e 2 
+ p—-Be ' 1 


(5.87) 


Inclination Angle Composite Solution. 


I c = (I)° + (I)* - [(I ) 0 ] 1 = (I)° + (I ) 1 - [ CD i ]° 

Combining Eqs (5.27), (5.59), and (5.70) gives a composite 

solution for the orbital inclination angle. 
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Combining Eqs (5.31), (5.59), and (5.66) gives a composite 

solution for the longitude of the ascending node. 
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(5.89) 


Argument of Latitude at Epoch Composite Solution . 

n c = (n)° + (n ) 1 - [(n ) 0 ] 1 = (n)° + (n ) 1 - [(n) i ] 

Combining Eqs (5.28), (5.58), and (5.72) gives a composite 
solution for the argument of latitude at epoch. 
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For Eqs (5.86) - (5.90), the constants are given by 
Eqs (5.79) - (5.83). Hence, solutions to the five equations 


















of motion have been developed as functions of , 

C , and C . 
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Matching £ Order Solutions to the dq/dh Equation of Motion 
An inner expansion of the outer expansion solution is 
undertaken as well as an outer expansion of the inner 
expansion solution. Matching is performed by equating these 
expansions of expansion solutions. 

Inner Expansion of the Outer Expansion Solution . The 
outer expansion solution for the dq/dh equation of motion 
was derived to order e and given by Eq (5.34). 
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Noting that £ = h/e and h = e £ and rewriting the 
above equation in terms of the inner variable, £, gives 
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Expanding for small e gives the inner expansion of the outer 
expansion solution to order e. 














[(q) 0 ] 1 = [C./2 - C C,)] 
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Outer Expansion of the Inner Expansion Solution . The 
inner expansion solution for the dq/dh equation of motion 
was derived to order e and given by Eq (5.61). 


(q) * = pr^-Be ^cos a + K + eK 
C D 1 7 


Rewriting this equation in terms of the outer variable, 
h, and expanding for small e gives the outer expansion of 
the inner expansion solution for q. 


[(q) 1 ] 0 = K - eK ? 


(5.93) 


Matching . The outer expansion of the inner expansion 
solution is equated to the inner expansion of the outer 
solution expansion. 
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Subtracting Eq (5.79) from this expression and dividing by e 


gives 
























The composite solution to the dq/dh equation can be 


easily constructed: 

q c = (q)° + (q) 1 - (q) 1 ] = (q)° + (q) 1 - [(q)°] 
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where C , C , and C „ are given by the following relations 
1314 ' 16 & & 


derived in Appendix C: 
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where 















and where C and C are constants of integration. 

6 2 0 


Hence, a solution for q has been found as a function of two 


constants of integration, C and C „, and the constants C , 

& ’6 20 ’ 13 ’ 


' , and C , which are dependent upon C and C . 

14 16 r 1 2 


Summary 


For convenience, a summation of the derivation results 
using matched asymptotic expansions are presented. 


Solutions to the Non-Rotating Earth Equations of Motion , 


Treatment of atmospheric entry as a boundary layer 
problem allowed for the application of the method of matched 
asymptotic expansions. In this section, the rotating Earth 
equations of motion were asymptotically expanded using a 
small parameter that is a function of the planet’s inverse 
atmospheric scale height, radius, and rotation rate. Zero 
order solutions to these equations were derived. These 
solutions were of two types, inner solutions, valid close to 
the Earth, and outer solutions, valid far from the Earth. 

The domains of these two sets of solutions overlapped. 

Hence, direct matchings of the expansions of the expansion 
solutions were accomplished. This produced equations that 
related the inner solution constants, K., to the outer 

solution constants, . Composite solutions, valid for the 

entire atmospheric entry domain of non-dimensional altitude, 


117 

























»'*■*■* * A. w » w *." A- 




were constructed using these relations. These composite 
solutions are solutions to the non-rotating Earth equations 
of motion for three-dimensional, lifting atmospheric entry. 


The outer solutions were given by 


I = C 
0 s 


= C „ 
0 4 


u (1+h) = C 
o v ' 1 


^ - V 1+h)2 


ri - c /(i+h) ■ 

a — - cos - 7 — + C 

0 . (1 - C 2 C )’ / 2 . 3 

v 1 2 ' 


(5.27) 


(5.28) 


(5.29) 


(5.30) 


u o = 1 + cos (a 0 - 0 3 )[l - c/cj 


2 n 1 1 /2 


(5.31) 


The inner solutions were given by 
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Matching the expansions of these solution sets produced the 

following relationships between the outer solution 

constants, C., and the inner solution constants, K.. 
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(5. 

The equations given above are for the three-dimensional 
atmospheric entry of a non-rotating Earth. These solutions 
were derived from zero order terms in the expansions of the 
rotating Earth equations of motion. 

The composite solutions for the three-dimensional 
atmospheric entry of a non-rotating Earth were given by 
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q c = Cj/[2(i+h) - c c 2 (i + h) 2 ] 


C D 


(5.87 


1° = cos 1 cosK cos tana 
4 


(5.88 


• log tan 4 + 2 


1 fCj , , 2 

^ + icos 1 Tr^Be ' E, cosa + K 

4 z l^D 1 


a = -cos 


1 - C/(l+h) 

.(1 - c 2 c 2 ) l/2 


+ sin sinK /cos cosK cos tana'log tan -7 
4 4 4 


1 -1 [ C L r -h/e 2 

+ „cos pr~Be 

2 l°D 


cosa + K 


+ cos 


1 - C 

_ 1 _ 

(1 - C 2 C )’ / 2 
v 1 2 J 


(5.89 


n c = K - cos 1 
5 


cos sin 1 fsinK /sin(I C )1 /cosK (5.90 

t 4 J 4 


Solution to the dq/dh Rotating Earth Equation of Motion 


To accurately model Earth atmospheric entry, a solution 


to the dq/dh equation that included rotating Earth effects 
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was required. This solution was derived in the form of 


order e inner and outer solutions to the asymptotically 


expanded dq/dh equation of motion for a rotating Earth. 


Matching was performed between these solutions and a 


composite solution to the dq/dh equation of motion for the 


three-dimensional atmospheric entry of a rotating Earth was 


constructed. 


The outer solution to the rotating Earth dq/dh equation 


of motion was given by 


o 

q = 


- c 2 (i + h) 


1 “I /2 


rc c c 

1 3 _ 1 4 _ , 14 

' <C 16 ) 2 ' C 1 6 


C exp -C h] 
2 0 r 1 6 J 


(5.34) 


The inner solution to the rotating Earth dq/dh equation of 


motion was given by 


q* = ^r^'Be ^cos o + K + eK 
°D 1 1 


(5.61) 


The composite solution to the dq/dh equation of motion is of 


order e and hence includes Earth rotational terms. The 


composite solution to the dq/dh equation of motion for the 


three-dimensional atmospheric entry of a rotating Earth was 


given by 
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where C „, C , and C were given by the following 
1314 16 0 & 

relations derived in Appendix C: 
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Traiectorv States of Validity for the Non-Rotatin 


Earth Equations of Motion 


Introduction 

In Section IV, the rotating planet terms in each of the 
five equations of motion for Earth atmospheric entry were 
examined. The dq/dh equation of motion for a non-rotating 
Earth was found to be never valid for the investigated 
ranges of inclination angle and speed. Because of this, a 
solution was developed for the dq/dh equation of motion for 
the three-dimensional lifting atmospheric entry of a 
rotating, spherical Earth. This was accomplished in Section 
V by treating atmospheric entry as a boundary layer problem 
and applying the method of matched asymptotic expansions. 

In Section IV, trajectory states were found to exist 
where some of the non-rotating Earth equations of motion are 
valid for a rotating Earth. The d(l/dh, da/dh, and dl/dh 
non-rotating equations of motion were all found to be valid 
for the same entry trajectory states. Other, independent 
trajectory states were found to exist where the du/dh non¬ 
rotating equation of motion is valid for a rotating Earth. 

In this section, these trajectory states are examined in 
more detail for the du/dh and dl/dh equations of motion. 
Plots of the Rotate term solutions to these equations are 
generated for a large range of realistic values of u, h, a, 
7 , and I. Trends in these solutions are discussed and the 


















non-existence of overlapping trajectory state solutions 
between the du/dh and dl/dh Rotate equations is verified, 
few methods to estimate solutions to the du/dh and dl/dh 
Rotate equations are also presented. 

Solutions to the du/dh Rotate Equation . In Section IV 
it was found that the rotating planet terms in the du/dh 
equation of motion are 


I Rotate 


l 1 /2 

= - 4 e 2 u(l+h)| cosl 


- 2 e 2 ( 1 +h) 


. 2 T . 

2 cosasin lsinq 
tan^ 


( 6 . 1 ) 


Setting this equation equal to zero and solving for flight 
path angle gives solutions for trajectory states where the 
rotating Earth terms in the du/dh equation of motion are 


7 = tan 

Rotate = 0 


2 

-l -Ccosqsinasin I 
2cosI 


( 6 . 2 ) 


where C = 1 

u 2 J 


(6.3) 


Solutions to the dfi/dh. da/dh, and dl/dh 


Rotate Equations . In Section IV, the solutions for the 
Rotate term expressions were found to be identical for the 
da/dh, dfl/dh, and dl/dh equations of motion. This result 
simplifies the search for trajectory states where the non¬ 
rotating Earth equations of motion are valid for a rotating 


Earth. Instead of detailed examination of the solutions to 













all three of the da/dh, dn/dh, and dl/dh Rotate term 
expressions, examination of only one of them is required. 
In this section the solutions to the dl/dh Rotate term 
expression are examined. From Eq (4.13) 


I Rotate 


= 2 fe iim 
2 u 


cosQsinl 
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(cosatan^ - sina) 
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2 utan^ 


(6.4 


Setting this equation equal to zero and solving for flight 
path angle gives solutions for trajectory states where the 
rotating Earth terms in the dl/dh equation of motion are 


. -l sina + Ccoslsina 

7 = tan - o - 

Rotate = 0 Zcosa 


(6.5; 


Non-Existence of Overlapping Solutions to the du/dh and 


dl/dh Rotate Equations 

To verify the non-existence of overlapping solutions to 
the du/dh Rotate and dl/dh Rotate equations, Eqs (6.2) and 
(6.5) are equated. 


-1 -Ccosasinasin I 


-1 sina + Ccoslsina 



















A trivial solution to this expression is given by sina = 0 . 
The non-trivial solutions to Eq ( 6 . 6 ) are found by dividing 
by sina and then solving for cosa. 


f-cosl 2 X 1 1 / • T 

cosa = I — - cos II /sinl 

„ h „ e c = 


(6.7) 


For real solutions to exist for Eq (6.7), the following 
condition must be true: 


cosl 

C 


+ 


cos 2 I < 0 


( 6 . 8 ) 


Since values of non-dimensional altitude and speed ratio 
are never negative, C always has a positive value. The term 
cosl is never negative for prograde entry trajectories. 
Hence, the condition given by Eq ( 6 . 8 ) is never satisfied 
and overlapping solutions to the du/dh and dl/dh Rotate 
equations do not exist. 


Graphical Trajectory State Examination 


















combinations of u, h, a, and I into the Rotate solution 
equations, given by Eqs (6.2) and (6.5), for the du/dh and 
dl/dh equations of motion. These solutions give trajectory 
states for which the non-rotating equations of motion are 
perfectly valid for rotating Earth entry analysis. This 
validity is dependent upon the physical characteristics of 
the planet under consideration, namely e ? , and the type of 

trajectory being analyzed; validity is not directly 
dependent on any of the vehicle’s physical characteristics. 
Determination of the validity of the non-rotating equations 
of motion for planets other than Earth could be undertaken 
in a similar manner as the analysis given in this section. 

A computer program was developed to solve Eqs (6.2) and 
(6.5) for a large range of u, h, a, and I. Plots of 
families of curves of 7 versus I are presented for various 
ranges of a and for realistic combinations of u and h. 
Appendix E presents three example Earth atmospheric entry 
trajectories to give an indication of what the dimensionless 
variables u and h are in relation to more conventional entry 
trajectory parameters such as Mach, altitude, and velocity. 
Most current and planned lifting entry vehicles have 
trajectories which fall within the ranges of u and h given 
by these examples. The realistic combinations of u and h 
discussed in the following pages are derived from liberal 
estimation of their ranges in these example trajectories. 
Values of u selected for investigation ranged from circular 












orbital velocities to low supersonic speeds where terminal 
maneuvers, such as landing approaches, are usually 
initiated. Realistic values of non-dimensional altitude 
corresponding to these values of speed ratio were used in 
this analysis. Values of argument of latitude at epoch and 
orbital inclination angle had ranges of 0.0 <_£»<_ 360 
degrees and 0.5 I £ 75 degrees, where most atmospheric 
entry occurs. Because of the large number of plots needed 
to display data trends, some of the plots which are 
discussed in this section are contained in Appendix F. In 
all the plots referred to in this section, a, the argument 
of latitude at epoch, is referred to as "Alpha". All values 
of 7 , a, and I given in these plots are in units of degrees. 

Trajectory States of Validity for the du/dh Non-Rotating 

Earth Equation of Motion . Trajectory states are 
presented below where the non-rotating Earth du/dh equation 
of motion is valid. Trends in these solutions to Eq (6.2) 
are discussed. 

IZR ical Solution Observations . Figure 13 is a 
typical plot of solutions to Eq (6.2) of flight path angle 
versus orbital inclination angle. For constant, non-zero 
values of argument of latitude at epoch, 7 decreases as I 
increases. Values of flight path angle also decrease in 


this plot as a increases. In addition, as a increases for a 
constant I, the change in 7 decreases, especially for large 
inclination angles. 




























13. 
















Effect of h Variation . Holding all other variables 


constant, the effect of change in non-dimensional altitude 
on these solutions is negligible. Figures FI and F2 in 
Appendix F present results for values of h higher and lower 
than in Figure 13 for the same values of u and a. 

Comparison of Figures FI and F2 indicates that only 
extremely small changes in flight path angle occur for this 
large change in h. This ineffectiveness of h variation on 7 
can also be intuitively seen from examination of Eq (6.3). 
Note that for all realistic values for non-dimensional 
altitude, 1.0 is much larger than h in the expression 

a * h ) 3 / 2 

Effect of u Variation . Figures F3 and F4 present 
solution results for values of speed ratio higher and lower 
than in Figure 13 for the same range of values of argument 
of latitude at epoch. Comparison of Figures 13, and Figures 
FI - F4 demonstrates that large changes in flight path angle 
occur for this large range of u. As u decreases, the 
magnitude of 7 increases for curves of constant a 

Effect of a Variation . The trend of the effect of 
variation in argument of latitude at epoch can be seen by 
review of F.gures 13 and Figures FI - F4. For the range of 
values of 0 in these plots, 7 and the change in 7 decrease 
as a increases along lines of constant inclination angle. 
These trends are more prominent for large values of I. The 
a terms in Eq (6.2) are 
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- cos(a)sin(a) 


Examination of this term and the properties of the inverse 


tangent function in Eq (6.2) indicates the solutions given 


in Figure 13 will be symmetric about a = 0 and 180 degrees 


and these symmetric solutions will repeat about 


±90 - a (degrees) for ± a 


Figures 13, F5, F 6 , and F7 demonstrate these 


observations. 



Trajectory States in Three Dimensions. Figures 14 and 


15 present three-dimensional plots of the trajectory states 


of validity for the du/dh non-rotating Earth equation of 


motion. Solutions to Eq (6.2) are plotted here for flight 


path angle versus orbital inclination angle and argument of 


latitude at epoch. The solution surfaces in these figures 


are constructed of lines of constant a and lines of constant 


flight path angle. These contours of constant 7 are drawn 



for changes of 7 of 0.165 degrees in Figure 14 and changes 


of 7 of 1.04 degrees in Figure 15. In both of these plots, 


values of flight path angle are near zero for small 


inclination angles and for very small values of argument of 


latitude at epoch. Although Figures 14 and 16 present 


solutions to Eq (6.2) for vastly different values of speed 


ratio, the two plots look almost identical; only the scale 


of flight path angle values differs. 


N 


1 




t 

I 


* 

i 


I 


I 












































































' States of Validity for the dl/dh Non-Rotatin 


Earth Equation of Motion . Trajectory states are given below 
where the non-rotating Earth dl/dh equation of motion is 
valid. Trends in these solutions to Eq (6.5) are discussed. 

Typical Solution Observations . Figure 16 is a 
typical plot of solutions of flight path angle versus 
orbital inclination angle. For constant non-zero values of 
argument of latitude at epoch, 7 increases as I increases. 
Values of flight path angle decrease in this plot as values 
of a increase. 


Effect of h Variation . Holding all other variables 
constant, the effect of change in non-dimensional altitude 
on these solutions is negligable. Figures F 8 and F9 present 
solution results for values of h higher and lower than in 
Figure 16, and for the same values of u and a. Comparison 
of Figures F 8 and F9 indicates that only extremely small 
changes in 7 occur for this large change in h. As in the 
case for the du/dh Rotate solutions, this ineffectiveness of 
h variation on 7 can also be intuitively seen from 
examination of the expression for C in Eq (6.5). 

Effect of u Variation . Figures F10 and Fll present 
solutions for values of speed ratio higher and lower than in 
Figure 16 for the same range of values of argument of 
latitude at epoch. Comparison of Figures 16, F10, and Fll 
indicates that large changes in 7 occur for this large range 
of u. As u decreases (along with the ineffectual non- 
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dimensional altitude), values of flight path angle increase 


for curves of constant a. 


Effect of a Variation. The trend in the effect of 


variations of argument of latitude at epoch can be seen by 


review of Figures 16 and Figures F 8 - Fll. For the range of 


a in these plots, 7 decreases as a increases for constant I. 


This change in flight path angle is more prominent for small 


values of speed ratio. Change in the range of a can also be 


easily seen analytically. The a terms in Eq (6.5) can be 


written simply as tan(a) . Examination of this term and 


the properties of the inverse tangent function in Eq (6.5) 


indicates the solutions given in Figure 16 will be symmetric 


about a = 0 and 180 degrees, and that these symmetric 


solutions will repeat about 


q * 180 (degrees) for * a 


Figures 16, F12, F13, and F14 demonstrate these 


observations. 


Trajectory States in Three Dimensions . Figures 17, 18, 


and 19 present three-dimensional plots of the trajectory 


states of validity for the dl/dh non-rotating Earth equation 


of motion. Solutions to Eq (6.5) are plotted here for 


flight path angle versus orbital inclination angle and for 


argument of latitude at epoch. The solution surfaces in 


these figures are constructed of lines of constant flight 


path angle and lines of constant argument of latitude at 


epoch. Figures 17, 18, and 19 display plots generated for 
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Figure 18. Surface of Trajectory State Solutions of Validity for the 
Rotating Earth, dl/dh Equation of Motion (u = .300 , h = 
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successively decreasing values of speed ratio. Although 
these values of u differ greatly, the solution surfaces are 
very similar in appearance. 


Solution Estimation Methods 

Some effort was spent in developing the computer program 
to solve Eqs (6.2) and (6.5) for large ranges of values of 
speed ratio, non-dimensional altitude, orbital inclination 
angle, flight path angle, and argument of latitude at epoch. 
To rapidly assess the validity of the non-rotating Earth 
atmospheric entry equations of motion for a particular 
trajectory state, there is an advantage to having a quick 
estimation algorithm. On the following pages, methods are 
presented to estimate some solutions to the du/dh and dl/dh 
Rotate expressions given by Eqs (6.2) and (6.5), 
respectively. 

Estimates to the du/dh Rotate Solutions . Figure 20 
presents a plot of trajectory state solutions to Eq ( 6 . 2 ), 
the du/dh Rotate equation. The solution surface in this 
figure is onstructed of lines of constant flight path angle 
and lines of constant argument of latitude at epoch. Two 
separate estimations of solutions to Eq (6.2) can be made 
with plots such as Figure 20. The first is the 
approximation that values of flight path angle in the 
triangular plane (ABC) of trajectory states have some 
constant, near-zero value. In general, this is a good 
approximation for 7 , especially for high values of u. The 




































ABC plane of "constant" 7 is bounded by 


0.0 <_ a <_ 45.0 degrees 
0.5 < I 75.0 degrees 

Line AB: a = -(0.60403)1 - (45.302) degrees 

Therefore, the conditions for assuming flight path angle has 
some near-zero value are given by 

(.02207)a + (.01333)1 < 1.00 (6.9) 

and by the ranges of inclination angle and argument of 
latitude at epoch given above. 

The maximum error in assuming 7 = 0.0 degrees in the 
ABC plane in Figure 20 is approximately equal to the value 
of 7 corresponding to the second flight path angle contour 
from the top of the plot. This maximum error in flight path 
angle is about 7 = -0.56 degrees for this plot. 

Application of this estimate to the corresponding ABC planes 
in Figures 14 and 15 give maximum errors in flight path 
angle of about 0.30 and 1.5 degrees, respectively. 

An estimate of a non-zero value for flight path angle in 
the ABC plane can reduce these maximum errors. Assuming the 
value of 7 in this plane is equal to half the 7 value of the 
second constant 7 contour in Figure 20 produces an estimated 
flight path angle of -0.28 degrees with maximum error of 
0.28 degrees. Likewise, estimates of 7 = -0.15 degrees 
and 7 = -0.75 degrees in the ABC planes in Figures 14 and 
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15 give maximum errors in flight path angle of about 0.15 
degrees and 0.75 degrees, respectively. 

Another approximation to the solution of Eq (6.2) is 
given by the simple assumption that flight path angle has 
near-zero values for small values of inclination angle. 
Maximum error in an estimated flight path angle is dependent 
upon the chosen value of inclination angle for this 
approximation. For example, assuming 7 = 0.0 degrees for 
I 15.4 degrees invokes a maximum error in flight path 
angle of about 0.28 degrees in Figure 20. Assuming 7 = 0.0 
degrees for I 29.2 degrees produces a maximum error in 
flight path angle of about 0.56 degrees in Figure 20. 

An estimate of a non-zero value for flight path angle in 
the ABC plane can reduce these maximum errors for the same 
range of inclination angle. Assuming 7 = -0.14 degrees 
for I < 15.4 degrees produces a maximum error in 7 of 
about 0.14 degrees in Figure 20. Estimation of a flight 
path angle value of -0.28 degrees for I 29.2 degrees 
gives a maximum error of 0.28 degrees. 

The accuracy of both of these estimations for solutions 
to the du/dh Rotate equation is heavily dependent on the 
value of speed ratio for a given trajectory state and the 
estimate of flight path angle. 

Planar Estimation to dl/dh Rotate Solutions . A planar 
estimation method was developed to obtain rough estimates of 
solutions to Eq (6.5). For given values of speed ratio, the 
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following relations give an estimate of the planar surface 
of solutions that appears in Figures 17 - 19. 

- I ' 1 - V (610) 

where 

c = -89.40 degrees 
I = 52.65 degrees 
a = -21.00 degrees 

and values of a, b, and c are polynomial functions of u, 
given by 

a - (106.701) - (271.537)u + (1209.57)u 2 

- (1840.30)u 3 + (895.29)u 4 

b = (5.657) - (90.002)u + (401.76)u 2 

- (611.66)u 3 + (297.66)u 4 

7 o = (-26.104) + (81.784)u - (365.55)u 2 
+ (556.77)u 3 - (271.00)u 4 

Examples of planar estimates of solutions to Eq (6.5), 
generated by the above method, are presented in Figures 21 - 
23. These estimated solutions correspond to the numerically 
generated solutions presented in Figures 17 - 19. A more 
exact method to estimate solutions to Eq (6.5) is given by a 
three-dimensional curvefit algorithm described below. 
















Figure 21. Estimation of Trajectory State Solutions of Validity for the Non 
Rotating Earth, dl/dh Equation of Motion (u = .860 , h = .0129) 
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Solutions. 


Using polynomial least squares curve fitting 


routines, a curvefit to the solutions of the dl/dh Rotate 
equation was developed for ranges of inclination angle and 
argument of latitude at epoch of 0.5 I < 75.0 degrees 
and -45.0 a < 0.0 degrees. This curvefit is three- 
dimensional in the sense that it predicts solutions for 7 
from Eq (6.5) for various I, a, and u. Since changes in h 
were shown to produce negligable changes in 7 , h was not 
included in the fit. The resulting algorithm produces 
accurate solution estimates. Figure 24 presents numerically 
generated values of 7 versus I, and Figure 25 presents 
corresponding curvefit generated estimations of 7 versus I 
for the same values of a, u, and h. 
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Figure 25. Estimation of Trajectory State Solutions of Validity for the Non 
Rotating Earth, dl/dh Equation of Motion (u = .020 , h = .0052) 



















VII. Conclusions and Recommendations 


Conclusions 

This effort investigated the validity of the non-rotating 
planet assumption for three-dimensional Earth atmospheric 
entry. This study was limited to entry at orbital 
inclination angles between 0.5 and 75.0 degrees, where most 
Earth atmospheric entry occurs, and vehicle speeds ranging 
from circular orbital velocities to low supersonic speeds, 
where terminal maneuvers (such as landing approaches) are 
usually initiated. Constant lift-to-drag ratio and 
ballistic coefficient were assumed along the vehicle’s 
flight path. Validity results are coordinate dependent 
since singu1arities exist in the equations of motion. On 
the basis of this investigation, the following conclusions 
are made: 

1. As a set, the five non-rotating planet equations of 

motion (Section IV) are invalid for Earth atmospheric 
entry. Hence, the non-rotating Earth assumption, 
common in analytical entry analyses, produces 
incorrect entry trajectory results. Solutions to the 
non-rotating Earth equations of motion were derived 
from expansions of the rotating Earth equations of 
motion using the method of matched asymptotic 
expansions. These solutions are given by Eqs (5.86) 

- (5.90). 
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The dq/dh equation of motion for a non-rotating Earth 
is never valid for the ranges of inclination angle 
and vehicle speeds mentioned above. Hence, a 
solution, Eq (5.95), was developed for the dq/dh 
equation of motion for three-dimensional lifting 
atmospheric entry of a rotating spherical Earth. 

This solution is valid for the entire non-dimensional 
altitude domain and it accounts for the Coriolis 
acceleration on the flight vehicle. The method of 
matched asymptotic expansions was employed in this 
derivation. 

A variety of realistic Earth atmospheric entry 
trajectory states exist where some of the non¬ 
rotating equations of motion are valid for a rotating 
Earth. The dn/dh, da/dh, and dl/dh non-rotating 
equations of motion are valid for the same entry 
trajectory states for a rotating Earth. Other 
trajectory states exist where the du/dh non-rotating 
Earth equation of motion is valid for a rotating 
Earth. These two sets of trajectory states do not 
overlap for the investigated ranges of entry 
inclination angle and velocity. 

A number of trajectory states are presented where the 
du/dh non-rotating equation of motion and the dfi/dh, 
da/dh, and dl/dh non-rotating equations of motion are 
valid. Holding speed ratio constant, realistic 
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variations of non-dimensional altitude causes 
negligible changes in these states of validity. The 
magnitudes of the flight path angles where these non- 
rotating Earth equations of motion are valid increase 
as speed ratio values decrease. Trajectory states of 
validity are symmetrical about various values of 
argument of latitude at epoch for the non-rotating 
Earth dl/dh and du/dh equations of motion. Two 
methods to estimate flight path angle values where 
these equations of motion are valid are also 
presented. 

Recommen dations 

Based on the observations of the assumptions and of the 
findings of this investigation, the following 
recommendations for further study are proposed: 

1 . Investigation of the validity of the non-rotating 
planet assumption for three-dimensional Earth 
atmospheric entry was limited to orbital inclination 
angles between 0.5 and 75.0 degrees. Further 
investigation should be undertaken for other ranges 
of inclination. Near-polar and negative inclination 
angles would be of particular interest. 

2. Further study of the trajectory states where the 


du/dh and the dfi/dh, da/dh, and dl/dh non-rotating 
Earth equations of motion are valid is warranted. An 
















attempt should be made to construct entire 
trajectories from these trajectory states. If 
successful, a reduction in the complexity of the 
equations of motion could be made for certain 
trajectories and possibly even specific trajectory 
classes. This reduction would allow the non-rotating 
planet equations of motion to be used in place of the 
more complex rotating planet equations of motion, 
causing past analytical studies that utilized the 
non-rotating planet assumption to be valid for these 
trajectories. Current trajectory optimization 
computer programs consume many valuable hours of 
computer run-time, iterativly integrating the 
rotating planet equations of motion. Optimization 
computer programs such as the Air Force Flight 
Dynamics Laboratory’s ENTRAN (ENtry TRajectory 
ANalysis) and OTIS (Optimal Trajectories via Implicit 
Simulation) codes, could greatly benefit from a 
reduction in the complexity of the rotating planet 
equations of motion. 

The method of matched asymptotic expansions has 
proven to be a powerful tool in the development of 
solutions to boundary layer problems such as 
planetary atmospheric entry. Further application of 
this method should be undertaken to obtain a full set 
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of solutions to the equations of motion for three- 
dimensional lifting atmospheric entry for a rotating 
Earth. This activity would entail finding solutions 
to the complex differential equations given by Eqs 
(5.20) - (5.24) and Eqs (5.46) - (5.52), and then 
matching and forming composite solutions from these 


results. 














Appendix A 


Derivation of Equations Relating 6, , j) and a, 0, I 

Using Spherical Trigonometric Relationships 


This appendix presents the basic derivation of equations 
relating 6 , <p, ifi and a, fi, I using trigonometric 
relationships. In this paper, the planetary model is 
spherical and hence the spherical sine and cosine laws can 
be successfully applied. 

Figure A1 depicts a reference spherical triangle. For 
clarity, the angles between the curved line segments are 
called interior angles whereas the angles formed between the 


unit vectors, r., r. , and r, , are termed exterior angles. 

1 J K 

Hence, A, B, and C are exterior angles and a, b, and c are 
interior angles in Figure A1. 

The sine law for spherical triangles states that the 
ratio of the sine of an interior angle to the sine of its 
opposing exterior angle is constant. 


sinCal sin(b) sin(c} , A , N 

sTn(A) = sin (B) = iK(C) (A1) 

The cosine law for a spherical triangle states that the 
cosine of an interior angle is the sum of two products. The 
first product is formed of the cosines of the other two 
interior angles; the second product is formed of the sines 
of the other two interior angles and the cosine of the 























Figure A1. Reference Spherical Triangle 

) 

s 
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opposing exterior angle. The cosine law is given by the 
following expression. 


cos (a) = cos (b) cos (c) + sin (b) sin (c) cos (A) 


(A. 2) 


Eqs (A.1) and (A.2) are used to relate the variables 
4> > 0, and I. From comparison of Figures 11 and 12 

(Section II) to Figure A1, the exterior and interior angl 
of the formed spherical triangle are 


e, 


i 

i 

i 


es 


Exterior Angles 

A = 0 - n 
B = 4 > 

C = a 


Interior Angles 
a = k/ 2 - f 
b = I 
c = k/2 
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where 


6 


longitude 

latitude 


<P = 

ip = heading angle 

a = argument of latitude at epoch 

0 = longitude of the ascending node 

I = inclination angle 

Application of the sine law to these variables gives the 
following: 


i 



! 

L' 

s 

Si 

V 

V 

I 

c 

* 




s_in^ _ sing _ sin(g - fl) 
sinl sin(7r/2) sin(7r/2 - ip) 

Since 


cosip = sin(^ - 

*> 

and 

sin(|) = 1 


sin(0 - 

sing =-^ — 

cos ip 

HI 



(A. 3) 

simp = sinlsing 



(A. 4) 


Application of the cosine law to the inclination angle 
gives the following. 

cos (I) = -cos (^) cos - r/>) + sin (^) sin cos (^) 

«■ «r 

Since sin^ = cos (g - ip) and cos (g) = 0 

cos (I) = cos(|t)cos(^) (A. 5) 

Another application of the cosine law gives 

cos(^) = -cos(I)cos(^ - ip) + sin(I)sin(^ - ^)cos(g) 
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which reduces to 


cos(I)sin(^) = sin (I) cos (y») cos (a) 


or 

Yet 


cosa = 

another application of 


tan^ 
tan I 


the cosine 


law gives 


(A. 6) 


cos(^ - i>) = -cos(|) cos(I) + sin (^) sin (I) cos (8 - [}) 

sin(^) = sin(I)cos(0 -(]) (A.7) 

Some other relationships can be found that are coupled 
to Eqs (A.3) - (A.7) and are useful for some occasions. 
Substituting Eq (A.5) into Eq (A.6) 

cosa = cos^cos(0 - n) (A.8) 

Substituting Eq (A.5) into Eq (A.3) gives 

sin < 9 - n > = IS! < AB > 


Rewriting Eq (A.6) and substituting in Eq (A.5) gives the 
expression 


sin (ip) 


cosasinl 


cos <p 

Substituting Eq (A.4) 


into this result 




The relations given by Eqs (A.3) - (A.10) are 
Section III to transform the equations of motion f 
of latitude, longitude, and heading angle to terms 
containing orbital inclination angle, longitude of 
ascending node, and argument of latitude at epoch. 


(A.10) 

used in 
rom terms 

the 
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sin(eb) - 


Therefore 


3,3 5,5 7,7 

e b , e b e b 

3 ! 5 ! ~ 7 ! + 


sin(eb) = eb + 0(e ; 

As an example, the following expansion of the cosine 

2 

function is derived to Q(e ). 

cos a + eb + e 2 cj = cos ( a ) cos |eb+e 2 cj - sin (a) sin Jeb + e 

= cos(a)^cos(eb)cos(e 2 c) - sin(eb)sin (e 2 c) J 

2 2 1 

- sin(a) sin(eb)cos(e c) + cos(eb)sin(e c) 

o 

= cos (a) - eb - sin(a) + 0(e ) 

The following relations were obtained through 
application of the Taylor series expansion: 

cos j^a + eb + 0(e 2 )] = cos(a) - eb’sin(a) + 0(e 2 ) 

sin a + eb + 0(e 2 )j = sin(a) + eb’cos(a) + 0(e 2 ) 

tan a + eb + 0(e 2 )l = tan(a) + ——- + 0(e 2 ) 

J cos 2 (a) 


tan[a 4 eb + 0(e 2 )]] = ^ 


tan (a) . 2 , >. 

' sin (a) 


4 - + 0 (£ S ) 


[ ein [* - eb * 0(e!) ]] = sdrno - £b ’— 1 Hfr 4 0(€l) 

' sm (a) 
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Appendix C 


Derivation of Solutions to Selected Outer Expansion 


Differential Equations 


Outer Expansion £ Terms 


The e terms of the outer expansions found in Section V 


are as follows: 


dh " (1+h) 


(C.l) 


% % 

(1+h)L u 


(C.2) 


(C.3) 


(C.4) 


dh (l+h)tan 7 


(C.5) 


Solutions to this set of differential equations are 


derived below. 


du/dh Equation Solution . Rearranging Eq (C.l) gives 


o -dh 


(1+h) 


or ln(u ) = -ln(l+h) + K 


Solution: 


u (1+h) = C 
0 v ' 1 


(C.6) 


_% V** "» "■ _% .*1 — *» % **» ~w j,”** — ^ fc > ■ I, * a •> | 
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dq/dh Equation Solution . From Eq (C.2) 


dq o % \ % ,1 

dh _ (1+h) u ~ 


Application of Bernoulli’s equation to this 
differential equation makes the solution derivation 
straightforward. Bernoulli’s equation is given by 

^ + P(x)y = Q(x)y n 


This form gives the solution (Beyer, 1984:315) 

r 

i/exp [ (l-n) Jp(x)dxj = (1-n) Q(x) [ (1-n) Jp(x)dxj dx + C 


where v = y n and n ^ 1 


lienee , + P ( h )% = Q C^) % 


where 


n = 3, v = q^ 


P(h) = iTThy ’ and 


ech) > h- = nW 

i o v y 


Solution: JyT+hy = ^- n (^ + ^) 

r 

q^ 2 exp -21n(l+h;j = ^ exp j^-21n (1+h) j dh + K 

1 - 


(l+h)l 


f 

• (1+h) _2 dh + 


















Integrating gives 


% 


V = § ( 1+h ) - c 2 (l + h ) 2 


(C. 


or 


% = [§ ( 1 +h) - C 2 ( 1 +h) 2 ] 


-1 / 2 


(C. 


dl 

_c 

dh 


dn c 

dh 


dl/dh 

Equation 

Solution. 

From Eq 

(C.3) 

= 0 ; 

I 

= C 





0 

5 




dn/dh 

Equation 

Solution. 

From Eq 

(C.4) 

= 0 ; 

n 

= C 





0 

A 




da/dh 

Equation 

Solution. 

From Eq 

(C.5) 


(C. 


(C. 


da. 


dh (l+h)tan 7 


By definition, q = cos'y . The expansions of q and C 0 S 7 
given by the following expressions: 

q Q + eq ] + 0 (e 2 ) = cos 7 q + £ 7 ^ 0 (e 2 )] 


= cos 7 0 - e 7j sin 7 0 + °(e ) 


By noting the order of e in this equation, it can be seen 
that 


7) 


8 ) 


9) 


10 ) 


are 


% = cos 7 0 


(C.ll) 









From Eqs (C.ll) and (C.12), an expression is found for the 
tangent of 7 . 


tan 7 = ——- 

'0 COS7 


(1 - cos ! 7o )’/ 2 (1-q, 2 ) 1 / 2 


tan To = [“T ” *] 


(C.13) 


lere was previously found to be 


= J-(l+h) - C 2 (l4h) 2 ] 


-1 / 2 


(C.14) 


Therefore 


k - c * (i+h)2 - i ]"' /2 -(ttet 


Making x = ( 1 +h) causes dx = dh and allows this integral 
to be rewritten as 


% = j[§. ' x - V 2 - 


-1 / 2 


The integral is now in the form of 


- dx 

°° = Jx(X ) ,/2 


(C•15) 


where X = a + bx + cx 2 , f - 4ac - b 2 . k = 4c/f. 
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i • 


f = 2 ^ C \ C 2 ~ 1) ' and k = “°1 C 2 /(C 1 °2 " 1} 

c 

1 


The solution to Eq (C.15) is 


a o = 


dh 


x (X) 


1 / 2 


= - (-a)' l/2 cos ! [ (hx + 2a)/{lxl (-q) 1/2 >] + < 


Hence 


cos 


-l 


2x/C - 2 


21x1 (C 2 -l/C 2 ) 1/2 


+ C. 


This solution can be expressed in either of the following 
two forms: 


a = - cos 


-l 


1 - C /(1+h) 


(1 - c 2 c 2 ) l/2 


+ c„ 


(C.ll 


or 


u = 1 + cos 


(a - C ) • f 1 - C 2 C 1 
v 0 3 J l 1 2 J 


1 / 2 


(c.r 


Outer Expansion e Terms 


The e 1 terms of the outer expansions found in Section 
IV are as follows: 


du -u 


l 1 


dh ” (1+h) 


] 1 /2 

3u_(1+h)I •cosl 


0 


0 


(C • 1 


dl 


dh 


= 2 


3(1+h) 


u 


1 / 2 


cosa sinl 


~ ~--(cosa tan 7 - sina ) (C.l 

tan7 v 0 '0 o' 
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dir + TTTKy = -4 cos(C 5 ).(3C 

This differential equation is in the form of 

* 

f 

^ 4 P (x) y = Q (x) 

This form gives the solution 
y' ex P [/p(x)dx] = Jq(x) [/p(x)dx]dx 4 C 
For the du^/dh differential equation 
y = u > X = h, 

P(h) = ~ (T+h) » and Q = - 4cos(C b ) * (3C ] ) 

exp [/(f;hy] = ex p[l n ( 1+h )] = (1+h) 

(l4h) U] = qJ (l4h)dh = |(l+h) 2 4 C 6 

Solution: 

C / 

u , = O^h) ~ 2(3C) 1/2 cos(C b )(l4h) 


dq/dh Equation Solution. 
























where u Q , u^ , and 1^ are given by the following equations: 


C /(1+h) 


(C.6) 


- 2(3C) 1/2 cos(C 5 )(l + h) 


(C.23) 


I = C 
o b 


(C.9) 


Substitution of Eqs (C.6), (C.9), and (C.23), into Eq (C.22) 


gives 


dq , V 

dh ~ Z C 


% q i 


2 (1+h) 


- 2q () (l+h)cosC & [|-] 


+ (3C ) ,/2 (l+h) 2 cosC 

C 1 1 b 


(C.27) 


iere q Q = §(l+h) - C ? (l+-h) 2 ] 


-1 / 2 


(C.8) 


Eq (C.8) can be substituted into Eq (C.27) to produce a 
complicated integral with h. However, for typical lifting 
earth atmospheric entry, altitudes of interest range from a 
minimum of 0 to a maximum of about 100 nautical miles. This 
is equivalent to a range of non-dimensional altitude of 
about 0.0 h < 0.029 . The maximum h encountered for a 

typical lifting entry trajectory for Earth is thus always 
less than the small parameter, e, allowing h to also be 
utilized as a small parameter in the outer expansion 
differential equations. The use of h as a small parameter 
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considerably simplifies the derivation of solutions to some 
of the first order e differential equations. 

Expanding the terms in Eq (C.27) for small h and 

neglecting the resultant terms of 0(h 2 ) greatly simplifies 

the analytical integration of this equation. One term that 

appears often in Eq (C.27) is q Q . Rewriting Eq (C.8) 


a = [p- + ^ - C - h 2 C - 2C hi 
o Lt t 2 2 2 J 


-1 / 2 


1 1 


A condition allowing the expansion of this equation for 
small h is given by 


§“(l+h) > C 2 (1+h) 2 


This condition can be easily seen to always hold true for 
any values of C^ , , and h. 

Expanding Eq (C.8) for small h using a Taylors series 
approximation gives 

% = [e; - C J ’ /! - h [e7 - c 2 ]'[§7 - 20 2 ] + 0(h2) < c - 28) 

Some other terms appearing in Eq (C.27) are likewise 
expanded: 


su( 1+h ) = 1 + 


o _ n n 


one 


r 

+ Oi 


(C. 29) 

















V 


= [§-- c 2 J - 3h [§- - c 2 ]/ IS- - C 2 J +0(h2) (c - 31) 


Defining two new constants, both dependent upon and C^ , 


simplifies Eqs (C.28) - (C.31). 


let C, J = [§-C 2 ] and C, , = [t- - C, ] 


With substitution, Eqs (C.28) - (C.31) become 


% = (C, ,) 


-1 / 2 


C h 

1 2 + 0 (h 2 ) 


(C ) 

v 1 l J 


l / 2 


(l + h) = (C ,)' l/2 + 


l 1 


c , < c ,, > 


3 /2 


+ 0(h 2 ) 


, 1 2C 1! h 2 
V = ^-- 0(h ) 


C M (C,,)’ 


V - < c ,,)' 3/2 - 7?TT72 + 0(h2) 


(C ) 
v 11 ' 


Two other terms in Eq (C.27) can also be expanded 


V 


(1+h) 2 = (C n ) 3/2 + 


(C ) 
v 1 1 J 


3 / 2 


3C 


2 - 


1 2 


1 1 




[TT = -q. + q.h + 0 (h 2 ) 


(C.32) 


(C.33) 


(C. 34) 


(C.35) 


(C.36) 


+ 0(h 2 ) (C. 37) 


(C. 38) 











Substitution of Eqs (C.33) - (C.38) into Eq (C.27) gives a 


solvable expression for dq /dh. 


<*q 


L = q, C, a + qC r h + C h + C 


dh ‘‘l 16 15” 14” 13 


(C. 39) 


where 


16 C (0 )'/ 2 C , C u 

1 v 11 J 


- 1 


(C. 40) 


2C 


c i6 “ 1 C C 


1 2 


1 1 1 


1 1 

C n + (C 1 ) l/2 


(C.41) 


3C C 
6 1 2 


14 (C,) 2 (C n ) 5 / 2 


+ 2cos(C ) 

D 


1-3 


1 2 


1 1 


c , 3 c ,, 3 


1 / 2 


(C. 42) 


C, 3 = 2cos(C 5 ) 


1 


3 




c c 1 


C C 

1 11 J 


l l l J 


T 1 /2 


3C 


C 2 C 3/2 

l i l 


(C.43) 


The differential equation given by Eq (C.39) is in the form 

+ P(x)y = Q(x) 

This form gives the following solution 
y ’ exp |/p(x)dxj = Jq(x)[/p(x) dxjdx + C 

For the dq^/dh differential equation y = q , x = h, 














With substitution 


exp[Jp(h)dh] = exp[-C ig h - §( c 15 ) h2 ] = ex p[ _C 16 h 4 °( h2 )] 


Therefore 


r r 

q i exp [~ C i6 h ] = c !3 ex p[-C 16 h]dh + C 4 h[exp(-C 6 h)]dh + C 


and the solution to the dq /dh equation is given by 


q i " C 


1 3 


1 4 


- h 


1 6 


(C )• 
V 1 6 ' 


1 4 


L C 1 6 J 


- C exp[-C hi 

20 K L )6 J 


(C. 


The constants C , C , and C are given by Eqs (C.43), 
1314 16 


(C.42), and (C.40), respectively. Because li terms are 


neglected, C is dropped from the solution. Hence, a 


l b 


solution for q has been found as a function of two 


constants of integration, C g and C^ Q , and the constants C 


C , and C , which are dependent upon C and C from the 
14 3 6 12 

relations given in Eq (C.32). 
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Appendix D 

Derivation of Solutions to Selected Inner Expansion 


Differential Equations 


Inner Expansion e u Terms 


The e° terms of the inner expansions found in Section V 
are as follows: 


o d 1 Lcosa 

0 sin 7 0 c d cos 7 0 


(D.l) 


^ 


Be * cosff 


(D.2) 




sinacosa 


cos 7 o sin 7o 


(D.3) 




sinasina 

_ o 

cos 7q sin 7o sinI Q 


(D.4) 


smasma 


tanl sin^ cosl 

0 '0 0 


(D. 5) 


Solutions to this set of differential equations are derived 


below. 


dq/d£ Equation Solution . From Eq (D.2) 


L n -f 

- ■ Be 'cosa 

C D 


This is a very simple integration by separation of 
variables. The constant of integration is K . 















(E 


a = rf — - Be Vosa + K 
^0 c D 1 


du/dh Equation Solution . From Eq (D.l) 


du 

37 ^- = - 2 u Be ^ 
d£ 0 


1 L.COSCT 


sin 7 0 c d cos 7 


0 J 


Previously it was seen that q^ = cos 7o 


T . , % d 1o . d 1o C L„ -f 

Therefore, ^7— = -sin).— and 3 = 7,—Be s 




0 d£ 


d S " C D 


The du/dh equation can be rewritten 


^“o - 


-2u Be 




0 sm) 


0 


-2u Be f 


0 


7) % 


d J±.d£_ 

d ? d ') 0 

d e 


fi.df, 
d « d % 


ln(u 0 ) = 


1 +1 




-2u Be 

0 Sin7 


i -1 


L. Be -£ cosa 
C D sm 7 


i +1 


-2u Be“« 

0 C D % J 


C L „ -5 

— Be '•cosa 
C D 


0 J 

-1 


d 7 r 


lq o 


After reducing, this equation takes the form 


ln(u 0 ) = 


f - C P 2d t + 

Cj^COSCT 


2 d 
% 


a = < y 2l ° + 21 n (7 ) + 

T) C^COSCT 'O 


Solution: 


u 0 = K 3% exp 


■ C D 2 7o 


Cj^cosct 


dl/d£ an d da/d£ Equations . From Eq (D.3) 


COSCT 

sin 7 0 


In ( K 3 > 

(D 


















a.L, o T . sinacosa 

—2. = _k. Be _ £-°- 

d i C D cos 7 0 sin 7 0 


From Eq (D.5) 


C T * sinasina 

- p—'Be £ • — — - ——j 

G~. tanl sinnr cos! 

D 0 '0 0 


By the chain rule 


dl dl d£ 

da, df da, 

0 s 0 


tanl 
_0 

tana. 


With separation of variables, this integration simply 
produces 


sina sinl = sinK 

0 0 4 


s i nK "* 

I = sin " 1 —--I 

0 sino. 


da/d£ and dn/d£ Equations . From Eq (D.5) 


^ sinasma 0 

tanl sinnr cosl 

0 '0 0 


From Eq (D.4) 


dn C T f. sinasina, 

o = _L. R -C- - --°_ 


d C ‘ % 


cos 7 q sin 7 o sinI Q 


Chain rule 


da Q da Q d£ 

dn d£ dn. 


= -cosl 


Substitution gives 


da, , r sinK ,1 r sinK , 2 ' 

j n = -cos sin [—:-1 = - 1 - [—t-1 

dn, lsina„J lsina A J 

0 . 0 . 0 


Rewriting this equation 














sina da 

0 0 


an = 


[sin 2 a 0 - 8in 2 Kj 


l / 2 


To solve this integral, let x = cosa and 


a = 1 - sin 2 K = cos 2 K , 
A A 


0 9 

so dx = - sina ’da and 1 - x = sin a 

0 0 0 


Therefore, the integral can be written in the form 


dx 


J [a 2 - x 2 ] 1 / 2 


= - cos 


■i r_x_i 

llal 


Solution: [} Q = - cos - K, 


or 


cosa = cosK cos(K - 0 ) 

U i o 0 


dl/d^ and d^y/d^ Equations . From Eq (D.3) 


dl C T .. cosa„ 

_0 q c • 0 

—— = *=—• Be s sina- 


d £ 


"D 


cos'v sin^ 


o 


From above 


d 7 n 


d £ 


C T f 

L. Be -C cosa_ 


Bin 7, 


o 


By the chain rule 


dI o dI o d * . cosa o 
d 7n d C d 7 n cos 7 n 


Since 


. -i r sinK <i 

a, = sin v--T~ 


then by substitution 









tana 


cos 7„ 


rsinK ,' 

cos sin 1 —:—?4- dl 

Lsinl J 0 

J L 0 J J 


This can be rewritten as 


tana 


COS7 


sin 2 K. 1 1/2 


• 2 T 
sm X. 


sinl 

_0 

. [sin 2 -1 + s 


hsin I 


2 T 11 /2 0 


To solve this integral, let x = cosl and 


2 i • 2 t/ 2 r/ 

a = 1 - sin K = cos K 
4 4 


so dx = - sinl dl and 
0 0 


1 - x 2 = sin 2 I„ 


Thus, the right side of the integral can be written in the 


and the left side of the integral is 

* j 

tan<T = Analog [tan [| + J-] ] + K g 

0 


Solution: 


cosI Q = cosK^cos[tana * log[tan + 
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Inner Expansion £ lerms 

The e 1 terms of the outer expansions found in Section V 
are as follows: 


C T - f- a , sina, 

p—-Be ^sina - - -- 

C D [ cos 7 0 sin 7o 


tana dl 
o . i 


'1 2 
COS 7 
»n 


r l • 2 
sm 


tanl d£ 


dOj s 

ar = cr Be tsl ™ %—? 

* U cos 


(D. 13 ) 


(D. 14 ) 


sma„ sma 

0 0 

-2- - 7j ~ 2 - 

cos 7 tanl sin 7 tanl 

*00 '00 


sma„ cosa„ 

- _ 0 _ + a _0_ 

"1 . 2 T costl sin 7 tanl 

cos 7o sin 7o sin I Q '0 '0 0 


o „ -£ cos 7 0 cosff 

= 2u » Be ■»,—— + eraser 

sm D '0 

•n 


! t C t 7 cosa 

L cosa L '1 


C r> • 2 

D sin 7 cos 7 


2u^ Be 


1 Lcosq 

sin 7 0 + C d cos 7 0 


(D. 15 ) 


dn tana dl 

_L _ _0 . 1 

d£ sinI Q d£ 


dn i C L „ -f . [ sina o 

ar = cr Bs sin<T t— 5 —rr 

s D cos 7^ sinl 

• O 1 


' 1 • 2 

sinl sin 7 sinl 

0 0 '0 0 


smasina, cosa, 

_ 1 _0_ + 7 ,_ 2 _ 

1 . . 2-r COS7, sin7 sinl 

cos7 0 sin7 o sm 1^ '0 '0 0 


(D. 16 ) 
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To solve the dq /dh equation is a trivial exercise; 


solutions to Eqs (D.13) - (D.16) are not required. 

dq/d^ Equation Solution . From Eq (D.17), q^ is simply 


q, = K ? 


(D.18) 











This appendix presents three example Earth atmospheric 
entry trajectories to give an indication of what the 
dimensionless variables u and h are in relation to more 
conventional variables such as Mach, altitude, flight path 
angle, and velocity. These trajectories were generated by 
AFWAL/FIMG, Air Force Flight Dynamics Laboratory personnel 
in the course of flight performance analysis support for 
current hypersonic entry vehicle studies. Modern trajectory 
analysis computer programs were utilized to produce this 
data. 

Example Trajectory 1 is the gliding entry trajectory 
for a relatively high lift-to-drag ratio vehicle (about 
3.0). Example Trajectory 2 is the gliding entry trajectory 
for a low lift-to-drag ratio vehicle (about 0.5). Example 
Trajectory 3 is a very steep gliding entry trajectory for 
the same low lift-to-drag ratio vehicle. This last 
trajectory is somewhat non-realistic due to the extremely 
high dynamic pressures (and hence severe heating) 
encountered. However, it is included to help give an 
indication of the extreme in u and h combinations. Many 
current and planned lifting entry vehicles have trajectories 
that fall within the ranges of u and h shown in these 
examples. 

















For the sample trajectories presented in Tables I, II, 
and III, the following variables are defined. 

y = Altitude (kft) 

7 = Flight path angle (deg) 

V = Velocity (kft/s) 
h = Non-dimensional altitude = r/r 


u = Speed ratio 
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f 


287.4 

262.3 

227.9 
189.6 
171 .0 

175.1 

164.1 

138.9 

129.4 
107.0 


Table II Example Traiector 


V fkft/s 


10.53 


h flO 


25 

.01 

13 

.73 

25 

.01 

12 

.53 

24 

.92 

10 

.89 

24 

.47 

9 

.06 

23 

.24 

8 . 

. 17 

21 

.93 

8 . 

.37 

tO 

O 

00 

t- 

7 

.84 

00 
1— < 

.29 

6 

.64 
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Table III Example Trajec 


f 


300.0 
256.5 
212.8 
169.3 


V fkft/s 


h no 


25 . 

.00 

14 . 

.34 

25 

.05 

12 

.26 

25 

.08 

10 

.17 

25 . 

.04 

8 . 

09 

22 

.48 

4 . 

.43 

18 

.59 

3 

.46 

10 

.76 

3 

.18 
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Figure F6. States of Validity for the Non-Rotating Barth 
du/dh Equation of Motion (u = .210 , h = .00' 
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Figure F9. States of Validity for the Non-Rotating Earth 
dl/dh Equation of Motion (u = .210 , h = .00 



















Figure F10. States of Validity for the 
dl/dh Equation of Motion 
























States of Validity for the Non-Rotating Earth 
dl/dh Equation of Motion (u = .010 , h = .00 
















Figure F12. States of Validity for the Non-Rotating Earth, 

dl/dh Equation of Motion (u = .210 , h = .0077) 
















Figure F13. States of Validity for the Non-Rotating Earth, 
dl/dh Equation of Motion (u = .210 , h = .007 
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